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Abstract 



The present paper is the companion of JjJ] in which we proposed a scheme that tries to 
derive the Quantum Field Theory (QFT) on Curved Spacetimes (CST) limit from background 
independent Quantum General Relativity (QGR). The constructions of Q make heavy use of 
the notion of semiclassical states for QGR. In the present paper, we employ the complexifier 
coherent states for QGR recently proposed by Thiemann and Winkler as semiclassical states, 
and thus fill the general formulas obtained in H with life. 



We demonstrate how one can, under some simplifying assumptions, explicitely compute 
*q | expectation values of the operators relevant for the gravity-matter Hamiltonians of Q in the 

complexifier coherent states. These expectation values give rise to effective matter Hamilto- 
nians on the background on which the gravitational coherent state is peaked and thus induce 
approximate notions of n— particle states and matter propagation on fluctuating spacetimes. 
We display the details for the scalar and the electromagnetic field. 
^ . The effective theories exhibit two types of corrections as compared to the the ordinary QFT 

on CST. The first is due to the quantum fluctuations of the gravitational field, the second arises 
from the fact that background independence forces both geometry and matter to propagate on 
a spacetime of the form R x 7 where 7 is a (random) graph. 

Finally we obtain explicit numerical predictions for non-standard dispersion relations for the 
scalar and the electromagnetic field. They should, however, not be taken too seriously, due to 
the many ambiguities in our scheme, the analysis of the physical significance of which has only 
begun. We show however, that one can classify these ambiguities at least in broad terms. 



1 Introduction 

Canonical, non-perturbative Quantum General Relativity (QGR) has by now reached the status of a 
serious candidate for a quantum theory of the gravitational field: First of all, the formulation of the 
theory is mathematically rigorous. Although there are no further inputs other than the fundamental 
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principles of four-dimensional, Lorentzian General Relativity and quantum theory, the theory pre- 
dicts that there is a built in fundamental discreteness at Planck scale distances and therefore an UV 
cut-off precisely due to its diffeomorphism invariance (background independence). Next, while most 
of the results have so far been obtained using the canonical operator language, also a path integral 
formulation ("spin foams") is currently constructed. Furthermore, as a first physical application, 
a rigorous, microscopical derivation of the Bekenstein-Hawking entropy - area law has been estab- 
lished. The reader interested in all the technical details of QGR and its present status is referred to 
the exhaustive review article and references therein, and to || for a less technical overview. For 
a comparison with other approaches to quantum gravity see |4|, || H . 



A topic that has recently attracted much attention is to explore the regime of QGR where the 
quantized gravitational field behaves "almost classical", i.e. approximately like a given classical 
solution to the field equations. Only if such a regime exists, one can really claim that QGR is 
a viable candidate theory for quantum gravity. Consequently, efforts have been made to identify 
so called semiclassical states in the Hilbert space of QGR, states that reproduce a given classical 
geometry in terms of their expectation values and in which the quantum mechanical fluctuations 
are small ]7|, || ||, 10, |llj. Also, it has been investigated how gravitons emerge as carriers of the 



gravitational interaction in the semiclassical regime of the theory [12, 13, 14]. The recent investigation 
of Madhavan and others [15|, 16], [17|, on the relation between the Fock representations used in 



conventional quantum field theories and the one in QGR further illuminate the relation between 
QGR and a perturbative treatment based on gravitons. 

In J]]] we developed and discussed a general scheme how one can define a theory of quantum 
matter coupled to quantum gravity in the setting of QGR and investigate its semiclassical limit. 
In the present paper we concretize the results of JTJ by employing a specific proposal H, jTUl ITT 



for semiclassical states for QGR. As the present paper relies on the general approach as well as on 
specific results of Q, it should be read together with the latter. Especially the discussion of the 
conceptual issues arising in the present context is much more completely covered in Also, it 
should be stressed that the cautionary remarks concerning our results made there apply even more 
to the present paper: The analysis of the semiclassical regime of QGR in general, as well as that of 
the coherent states ||, |l(J, [lTJ for QGR specifically has only begun recently, and so the main purpose 
of our work is exploratory. 

In the present paper, we roughly proceed in three steps: Firstly, we review the coherent state 
family introduced in || [TTj, |TT|, [TJ5] and fix the parameters in its definition in such a way that best 
semiclassical behavior is obtained for the observables relevant to our considerations. Then, under 
some simplifying assumptions, we compute the expectation values in the coherent states for the 
operators relevant for setting up the effective QFT for the matter fields according to J1J. Finally, we 
use the resulting effective theory to approximately compute the quantum gravity corrections to the 
dispersion relations for the scalar and the electromagnetic field. 

Let us consider these steps in more detail: 



In P, [10], 0, |TJ 



a promising family of semiclassical states have been constructed and analyzed. 
Each member of this family is labelled by a (random) graph 7 and a point meMin the gravitational 
phase space. Other states derived by the complexifier method |2(J could be used as well but for the 
exploratory purposes of this series of papers it is sufficient to stick to those simplest ones. 

Three scales enter the definition of the coherent states and are of considerable importance for their 
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semiclassical properties. These scales are the microscopic Planck scale £p, the mesoscopic graph scale 
e which represents the average length of an edge of 7 as measured by the three metric determined 
by m and a macroscopic curvature scale L which characterizes the scale at which matter (and thus 
geometry) vary. While £p,L are determined by the input m, the scale e is a priori a free parameter. 
We fix it by asking that a natural family of observables be well approximated by our coherent states 
which leads quite generically to a geometric mean type of behavior, concretely e oc l p L}~ a where 
< a < |. In contrast to the weave proposal |7| the graph scale is larger than the Planck scale 
due to the fact that we do not only approximate the three geometry but also the extrinsic curvature 
which forces the coherent state to depend on all possible spin representations of SU(2) and not only 
the defining (or any other single) one. 

The analysis of |T3 revealed that the coherent states proposed do not approximate well coordi- 



nate dependent observables like the holonomy or the electric flux operator. However, we discovered 
that operators which classically correspond to integrals of scalar densities of weight one are extremely 
well approximated. This class of observables contains Hamiltonian constraints and all spatially dif- 
feomorphism invariant quantities which suffice to separate the points of the diffeomorphism invariant 
phase space. The intuitive reason for this is the following point which has been stressed for years, 
among others, especially by Rovelli [21 , 22| : Matter can be located only where Geometry is excited! 



Classically this follows from Einstein's equations. In the quantum theory it is reflected by the fact 



that matter and geometry degrees of freedom are necessarily located on the same graph [^3| ^4 
Imagine now constructing a diffeomorphism invariant area operator Ar. In contrast to its companion 
Ar(5) well studied in the literature it does not depend on an externally prescribed coordinate surface, 
rather in order to model the measurement of the are of the desk table on which you are working 
right now one would construct a coherent state of the combined matter and geometry Hilbert space 
which is peaked on flat space and, say, on an electromagnetic field which is zero everywhere except 
for a region in the vicinity of the table. This way the dynamics automatically forces the surface to 
be adapted to the graph on which the coherent state depends. 

In a next step, we compute coherent state expectation values for the gravitational degrees of free- 
dom that appear in the matter-geometry Hamiltonians. This computation, although straightforward 
in principle, turns out to be quite tedious in practice. To keep the computational effort on a tolerable 
level and maintain some clarity of presentation, we simplify things by doing the calculation only for 
the Abelian (Idnii-Wigner) limit U(l) of SU(2) as gauge group. The computations done in ITOL ITT 



exemplify that this replacement does not change the results qualitatively, and therefore seems accept- 
able for the exploratory purposes of the present paper. A calculation in full generality should only 
be carried out after other issues have been settled, and will probably necessitate the use of computers. 



Ground-breaking work on the phenomenology of QGR has been done in [25. [26|. ^7], 28 1. In 
these works, corrections to the standard dispersion relations for matter fields due to QGR have been 
obtained. Since we are dealing with a theory for matter coupled to QGR in [Tj and the present work, 
it is an important question whether these results can be confirmed in the present setting. Therefore, 
as a final step, we formulate effective matter Hamiltonians on a graph based on the expectation 
values obtained before. The resulting theory is that of fields propagation on a random graph. It 
bears a remarkable similarity to models considered in lattice gauge theory |2T| [H]], and there is 



also a close analogy to the propagation of phonons in amorphous solids. As we have discussed at 
length in the resulting dynamics for the matter fields is very complicated, and analytic results 
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in the literature on lattice gauge theory and on amorphous solids are few. (To say the least. See 
however |32| for a beautiful numerical study of some two dimensional models from condensed matter 
physics.) Already a simplified one dimensional system (whose definition along with some results 
was sketched in |lj and will be covered more completely in |33[]) shows many of the complications 
(optical and acoustic branches, fuzziness of dispersion relations at high energies etc.) that are to 
be expected for the dynamics of fields propagating in a QGR background. Therefore, to compute 
dispersion relations for the models obtained, we have to rely on an approximation scheme denoted 
"graph averages", geared to the description of the dynamics in the limit where the energy of the 
fields is low (or, equivalently, their wavelength large). This approximation scheme leads to precise 
numerical values for all correction coefficients in the dispersion relations, once we have fixed a random 
process that generates our sample graph. The results we obtain are similar to those of p5| , p7|j in 
many respects, but differ in the scaling of the corrections. 

The validity of the approximation scheme we use has been discussed in [|I] but certainly merits 
future investigation. In any case, the resulting formulas can probably effectively handled in full 
generality only by a computer. 

Let us finish with a brief description of the contents of the sections to follow: 

The next section contains a short review of the construction of the complexifier coherent states 



10, 11 



In section |] we analyze the relation between the different scales that enter the definition of the 
coherent states, and their semiclassical properties. Relying on this analysis we fix the parameters of 
the coherent states for the rest of the paper. 

Section [| is the longest of the present paper. We show how expectation values in the coherent 
states can be computed and do the concrete calculations for the operators occurring in the Hamilto- 
nians for scalar, the electromagnetic and fermionic fields coupled to gravity. 

Section |^ deals with the computation of dispersion relations. We implement the procedure out- 
lined in |lj and compare our results to the ones in the literature. 

Finally, in section ||| we summarize what we have tried to do and what could be achieved with 
present technology. We conclude with a list of the open conceptual and technical questions that this 
work has left us with. 



2 Complexifier Coherent States 

The purpose of the present section is to review the construction and basic properties of the coherent 
states for QGR ||, [K], TT| . For an introduction to the formalism of QGR as a whole we refer the 
reader to 0, 3, or to the brief introduction in [|TJ. 

As already pointed out in the introduction, the task of constructing semiclassical states for 
QGR has received much attention J7|, [| [TO, (TTJ. Semiclassical states are states, so far in the 



kinematical Hilbert space of QGR, that approximate a specific classical geometry in the sense that 
expectation values of observables in such a state are close to the respective classical values and the 
quantum mechanical fluctuations are small. These requirements can certainly not be met for all 
possible observables, so the definition of a semiclassical state also involves specification of the class 
of observables that are well approximated. 
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In the present work, we will use the gauge theory coherent states (GCS for short) constructed in 
|| and subsequently analyzed in ]10|, [11], [HJ. These states are only one example of a large class of 
semiclassical states, called the complexifier coherent states. We refer to [EO] for an investigation of 



this class of states as well as a discussion of the relationship to ||, [15|, [L6|, [T7|, jlBf . 

The main mathematical tool used in the construction of the GCS is a generalization due to 
Hall [|34|, 3f| of the well known coherent states for the harmonic oscillator. The basic observation 



underlying this generalization is that the harmonic oscillator coherent states can be obtained as 
analytic continuation of the heat kernel on M. n : 



^ z (x)=e- tA 5 x ,(x 



x E M n , z G C n , 



the Laplacian A in the above formula playing the role of a complexifier. 



It has been shown in |34| that coherent states on a connected compact Lie group G can analogously 



be defined as analytic continuations of the heat kernel 

^ g (h)=e- tA -6[?\h) (2.1) 

to an element u of the complexification G c of G. These states have nice mathematical properties. 
Among other things, they are minimal uncertainty states for a certain pair of operators and they 
form an overcomplete set in the Hilbert space over G derived from the Haar measure. 

The case of this construction relevant for the definition of GCS is G = SU(2). Its complexification 
is given by SL(2,C) and can be parametrized as 



u 



exp \irjp 3 /2]h, p k £ 1R 3 , h E SU(2) (2.2) 



where ir^, k = 1, 2, 3 denote the Pauli matrices. 

A crucial question in view of applications to the construction of semiclassical states for QGR 
is whether the states ( |2.1| ) obey peakedness properties analogous to that of the harmonic oscillator 
coherent states. In |10] it was shown that this is indeed the case: For u given by p, h via the 
parametrization (|2.2j), the following holds: 



• %l) l u is exponentially (Gaussian) peaked with respect to the multiplication operator h on the 
group at the point h. The width of the peak is approximately given by \ft 

• ipu i s Gaussian peaked with respect to the invariant vector-fields at a point p/t in the associated 
momentum representation. The width of the peak is approximately given by l/\/t. 

For a more precise formulation of these statements we refer to |l(| . 



In QGR, the configuration degrees of freedom are represented by the holonomies along edges e of 
a graph 7 embedded in S. To use the coherent states on SU(2) for the construction of semiclassical 
states for QGR, momentum observables, that are associated to a graph in a similar way as the 
holonomies have to be defined. This was done in |3lfl . The construction can be summarized as 
follows (for the many details we refer the reader to the original work): To each graph 7 fix once and 
for all a dual 2-complex P 7 , i.e. roughly speaking a set of surfaces (S e ) eg £( 7 ) which intersect each 
other in common boundaries at most and such that the edge e of 7 intersects only S e and that this 
intersection is transversal. The surfaces S e shall be given an orientation according to the orientations 
of the edges e, i.e. the pairing between the orientation two form on S e with the tangent vector field 
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on e at the intersection point should be positive. Also to each point p lying in a surface S e fix an 
analytic path p(p) connecting the intersection point S e fl e with p and denote the part of e from e(0) 
to S e n e by e in . 

With the help of these structures, we can now define the quantity 



(2.3) 



where a is a length scale introduced to make Pj dimensionless and whose relation with tist = £ 2 P /a 2 . 
The key feature of this new variable is that 



-2 ^ , Jl 



{Pph e ,} = 'Ls^Tlhef, {Pi,Pj} = -^e^ijkPk (2.4) 

where e^fc are the structure constants of SU(2). Therefore, if h e is represented by the multiplication 
operator h on the cylindrical subspace corresponding to e, can be represented by the right invariant 
vector- field itX^ acting on the cylindrical subspace corresponding to e. 

Having the momentum variables p? at disposal, the construction of the GCS can now be finished. 
It needs three inputs: 

• A point (A(°\ E^) in the classical phase space that should be approximated. 

• A graph 7 and a corresponding dual polyhedronal decomposition P 7 of E, and the associated 
path system Il 7 . 

• The parameter t or, equivalently, the length scale a. 

For each edge e of the graph 7, one can now compute the holonomy in the classical connection 
and the classical quantities pf^ depending on A^°\E^ as expressed in (|2.3|) . The gauge coherent 
state for QGR is then defined as 

N 

(h 

71=1 

where e\, . . . represent the edges of the graph 7 and g e {Ao, E ) = exp(p^ e Tj / '2)he ■ 

The states thus defined inherit the peakedness properties of the coherent states (|2.1|) in an obvious 



way with respect to the elementary observables h ei , . . . , h eN and p^ 1 , . . . ,Pj N ■ For more complicated 
observables, a more detailed consideration has to be given. This is the topic of the next section. We 
will see that this analysis fixes the parameter t as well as the average edge length of the graph G, 
thus reducing the freedom in the construction of the GCS considerably. 

3 Observables and Scales 

In the previous section we saw that the complexifier coherent states VV7 wm ^ e use d in this 
paper (see [^0| for generalizations) depend on a point m G M. and a triple (7, P 7 , Il 7 ) where 7 is a 



graph, P 7 a polyhedronal decomposition of £ dual to 7 and II 7 is an associated path system. The 
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states 

^ 'in '"y ^1*6 linear combinations of spin network states over 7 (and all of its subgraphs) with 
coefficients which depend on m, P 7 , Il 7 . We are interested in the question which kind of operators O 
are approximated well by these states, that is, for which holds that expectation values are close to 
the classical value and for which the fluctuations are small. 

By construction, they approximate very well the holonomy operators h e and the electric flux 
operators Ej(S e ) where e runs through the set of edges of 7 and S e is a face in the polyhedronal 

decomposition dual to e. But how about more general operators such as h p ,Ar(S) where p is an 
arbitrary path and S an arbitrary surface? First of all, unless p is a composition of edges of 7 we 
have < il) mr f, h e ij) mr1 >= due to the orthogonality of spin- network states. Secondly, the expectation 
values of the area operator suffer from the "staircase problem" which says that unless S is 



composed of the S e then its expectation value will be off the correct value. 

The first reaction is: The states are not good, they must be improved. One such improvement 
could be by averaging over an ensemble of graphs || but as shown in this still does not improve 
the holonomy expectation values. Thus, one could think that one should construct semiclassical states 
of a completely different type, maybe going to a new representation of the canonical commutation 
relations Jl5|, [16], [lTj. However, this is not easy if the present formulation of QGR is to be kept as 



shown in PD |. It therfore seems that we are in trouble. 

There is a second possibility however: Maybe we are just trying to approximate the wrong 
observables? Notice that it is a physical input which observables should be well approximated, 
certainly we do not expect all classical quantities to be approximated well in the quantum theory. 
This is even true for simple finite dimensional systems such as the harmonic oscillator: The energy 
itself is well approximated but not its exponential. In our case, traces of holonomy operators and 
area operators are certainly natural candidates for operators to be well approximated because they 
are gauge invariant, suffice to separate the points of the gravitational phase space and are simple 
functions of the basic operators that the whole quantization is based on, namely holonomy and 
electric flux operators. Is it possible that there are observables which are better suited for our 
semiclassical considerations? 

A first hint of how such observables should look like comes from the observation that the volume 
operator Vol(i?) for a coordinate region R does not suffer from the staircase problem. A detailed 
analysis shows that this happens because the region R corresponds to a three-dimensional submanifold 
of E rather than one - or two dimensional ones. We therefore are led to the proposal that one should 
not look at holonomy and area operators but rather at quantities that classically come from three- 
dimensional integrals. There are classical observables of that kind that one can construct and which 
separate the points of the gauge invariant gravitational phase space as well: Let ui a be a one form, 
say of rapid decrease, and consider 

r E a E b 

Q(u) := / d z x^=J=uj a uj h (3.1) 
Jt. ydet(g) 

r B a B b 

M{u) := / d 3 x 3 3 f UgUb (3.2) 
Jt, y/aet(q) 

where EjEj =: det(q)q ab and = \^ abc F^ c and where F is the curvature of the connection A. Notice 
that both ( |3.1| ), ( |3.2| ) are of the type of operators that can be quantized with the methods of |23] in 
a background independent fashion since they are integrals of scalar densities. Moreover, they suffice 
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to separate the points of the gauge invariant phase space as one can see by suitably restricting the 
support of uj a and by the polarization identity for quadratic forms. 

The crucial fact about these quantities is now as follows: When we quantize them along the lines 
of p|J] they become diffeomorphism covariant, densely defined, closed operators on the kinematical 
QGR Hilbert space TLq of the following structure 

6{uj)T s = 22 ui(e)u(e')O v;e , e ,T s =: ^ O^u, v)T s (3.3) 

where T s is a spin-network state with underlying graph 7(s) and V(7), E{^) denote the sets of vertices 
and edges of a graph respectively and o;(e) = J u. The fact that an action only at vertices takes 
place in (|3.3|) is due to the appearance of the volume operator which enters the stage due to the factor 



of 1/ T^deU^) in ( |3.1| ) , ( |3.2| ) which is required by background independence and the requirement that 



only density one valued quantities can be quantized in a background independent way p4| . The 
operator O v . Bie i is a polynomial formed out of holonomy operators along the edges 0/7(5) and powers 
of the volume operator restricted to an arbitrarily small neighborhood of the vertex v. Now the 
coherent states are constructed precisely in such a way that the holonomy operators along the edges 
of 7(3) are well approximated and, as we will explicitly prove in this work, they also approximate 
very well the volume operator of [37|, |38| at least if the graph is six-valent, e.g. of cubic topology. 
(For other graph topologies the prefactor which enters the square roots that defines the volume 
operator, would presumably need to be adapted to the vertex valence, it should be larger (smaller) 
for valences smaller (larger) than six). 

Thus, due to the Ehrenfest properties proved in [|TTJ we conclude that at least for coherent states 
based on graphs with cubic topologies the operators (|3.1|), ( |3.2| ) are approximated well (with small 
fluctuations) provided the expectation values of ( |3.3|) define a Riemann sum approximation of the 
classical integrals ( [3.1|) , (|3~2l ). This is, however, the case by the very construction of such operators 



as outlined in [E3|. Thus, the mechanism responsible for the fact that no such problems as for the 
area and holonomy operators arise is due to the fact that for operators coming from volume integrals 
the elementary electric flux and holonomy operators involved are automatically those adapted to the 
graph in question. 

That only cubic graphs should give rise to the correct classical limit might be disturbing at first 
but it is on the other hand not too surprising: The volume operator at a given vertex v is a square 
root of an operator which in turn is a sum of basic operators, one for each unordered triple of distinct 



edges incident at v in |37| and in []38[ one considers only those triples which have linearly independent 
tangents at v. Each of these basic operators is a third order homogeneous polynomial in electric flux 
operators. With respect to our coherent states, each polynomial gives a contribution of the same 
order of magnitude. If n is the valence of the vertex of v then there will be altogether N(n) terms 
where N(n) = n[n — l][n — 2] for the operator in |J7| while for the operator of this number 



is smaller whenever there are triples of edges with co-planar tangents at v. The smallest valence 
for which the volume operator does not vanish is n = 3 in which case N(3) < 6. Since each term 
corresponds to the volume of the cell of the polyhedronal decomposition dual to 7, the factor 1/48 
dividing the sum over triples is too large. Now iV(4) < 24 is still too small while N(5) < 120 is 
already definitely too large for the volume operator of p7| . For the cubic topology we have, however, 
precisely N(6) = 48 for |38j because the only triples that contribute are formed by those spanning 
the eight octants defined by the coordinate system defined by the tangents of the six edges at v. For 
graphs of higher valence, unless there are sufficiently many coplanar triples, the [B3] volume operator 
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also over-counts the classical volume. Notice that none of these statements proves that one operator 
is proved over the other, it just means that our coherent states do not approximate both equally 
well. Only if one would know that our states are "the correct choice" , could one distinguish between 
the two kinds of volume operators on physical grounds. 

Intuitively, it is actually not too bad that only graphs of low valence should give rise to the 
correct classical limit. After all, one would not try to approximate a classical integral by Riemann 
sums in terms of graphs with vertices of arbitrarily high topology. Such graphs should describe 
quantum states without classical correspondence. It is also natural that cubic graphs are somehow 
distinguished because the classical integral is locally defined by a Cartesian coordinate system. 

Having convinced ourselves that the coherent states of the previous section actually do make 
sense at least for operators of the kind ( |3.1| ) and ( |3.2| ) we turn to the question how the scale e should 
be chosen. In order to quantize the classical integral 

Q(u) = [ d?xO ab u a u b (3.4) 



the procedure adopted in |24j was to define the operator on the spin-network basis. Thus, let 7 be 
a graph and v 1— > R v a partition of S where v runs through V(7). Let el be the coordinate volume 
of R v . Then we have 

0(u)= I ^O a W~^e30 ab (t;K(t;V f) (t;)=:0 7 M (3.5) 

veV{-y) jRv v 

where in the last step we have replaced the integral by a Riemann sum. The quantization of the 
term at v in the sum in ( |3.5| ) gives rise to the operator 

6 7 (u,v)= ^2 u(e)uj(e')d v - efi > (3.6) 

uede,de';e,e'e-E(7(s)) 

in (|3.3|) and by construction its expectation value in a coherent state ipy, m gives back 7 (a>, v)\ m to 
zeroth order in h. Thus, apart from quantum corrections for the expectation value, which we will 
call a normal ordering error, already the quantity 0(uj,m) — 7 (u; m ) is in general non-zero. This 



classical error error will decrease with e. With the Euler-MacLaurin error estimation methods of |L9 
one can prove an estimate of the form 

uj,m) (3.7) 

Lj 

where (3 > 2 and L is the average size of the the quantity [O ab u a ujb]" /[O ab u a uJb] where the double 
prime denotes second derivatives. Thus, L captures information about the gravitational curvature 
as well as the curvature of u. The size of (3 depends strongly on the randomness of the graph in 
question and also would change if one would average over graphs. 



More precisely, if we are interested in diffeomorphism invariant quantities (|3.1| ), (|3.2| ) such as the 
matter Hamiltonians that we wish to approximate in the following sections, then we should set, e.g., 
u a = 4>, a where is a scalar field or we should consider integrands of the form q ab E a E b / y det(g) 
where E a is the Maxwell electric field. To see what the matter and geometry scales involved are, 
consider the time - time component of the Einstein equations for electromagnetic waves with vector 
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potential A = ^e**'*!* - **). If q 2 is the electric charge, then the matter energy density is of the 
order A\k 2 jq 2 . If R denotes the curvature radius of the curvature tensor then we get from Einstein's 
equations R~ 2 ~ (£pA k) 2 /a where a = hq 2 is the Feinstruktur constant. Thus, if we introduce the 
wave length by k = 1/A then R~ 2 ps (10Ao£p) 2 X~ 2 <C A -2 at least for weak electromagnetic waves 
Aq <C 10 32 cm _1 . Thus, L should, for the applications of this paper to be thought of being very close 
to the matter wave length A and R is large, so that the geometry is almost flat. 

Let us now consider fluctuations. Since the quantities M(u), Q(oS) have different physical units, in 
order to compare their fluctuations we should compare their relative fluctuations which are dimension- 
free quantities. More precisely, we consider the expectation value in the coherent state if) 1)Tn of the 
relative deviation squared [0/0(m) — l] 2 between the operator O and its expected classical value 
0(m) which is a proper measure for the total deviation of the operator from the classical quantity 
due to 1) the fluctuation of the gravitational field and 2) its discrete nature which forces us to work 
with graphs rather than continuous integrals. If we denote by < . > m , 7 the expectation value in the 
coherent state ip mn and if there is no normal ordering error then we arrive at 

< - I] 2 >m^< [T^T - I] 2 >m. 7 - l ? ( 3 - 8 ) 

U{m) Uj{m) u{m) 

The second term in Q3.8D is of order (e/L) 2/3 as derived above. Now we see that for the quantity 
M(u>) the first term is divergent for flat data because M(u) = while for O = Q(u) there is no such 
problem. This is like comparing the relative fluctuations of x,p for the harmonic oscillator at the 
phase space point (x,p) = (0, 1) which of course makes little sense. To deal with this problem we 
chose the following strategy: We compare the relative fluctuations at generic points in phase space 
where we find a relation between the scale a of the coherent state and the scale L and then extend this 
relation to all points in M.. This strategy is certainly ad hoc but we do not see any other possibility 
at this point to fix the size of a by a more physical requirement. 

Accepting this we will consider non-flat data in which case generically L m R is closer to the 
curvature scale. If we assume that the operators 0^(uj,v) in ( |3.3|) are much weaker correlated for 
distinct v than for coinciding v (as it turns out to be the case) then we obtain 

< 6 2 > mn -(< 7 > m , 7 ) 2 w J^[< o 2 7 > m>7 -(< 6 V)1 > m , 7 ) 2 ] (3.9) 

V 

Restricted to 7, the operator O vn is a homogeneous polynomial of some rational power of the op- 
erators Pj ~ Ej(S e )/a 2 which are of order E e 2 /a 2 where E is some average value of and a is 
the coherent state scale introduced in the previous section. It is also a polynomial of some integral 
power of the operator h a — h~ l which is of order B e 2 where B is some average value of Bj and 
is approximately given by EqL~ 2 . As shown in ||19|| , the fluctuations for the respective vertices v 
is effectively given by exchanging JjV (m) by td 2 ltV (m)/[dP(S e )} 2 ~ t0 lsV / P(S e ) 2 for the electric 
fluctuations and by t9 2 /i 2 7 ^(m) /[dh a ] 2 ps t0 7>v /[h a — h^ 1 ] 2 for the magnetic ones where a is some 
loop incident at v. Inserting P(S e ) « E e 2 /a 2 and h a — h a m E e 2 /R 2 w E Q e 2 /L 2 . equating ( |3.9|) for 
O = M(u), O = Q{uj) respectively immediately leads to a w L. 

While the derivation of this result is maybe not entirely convincing, it is actually the only choice 
from a classical point of view: Since L is the only classical scale available and the complexifier 
generator C for our coherent states, from which the scale a derives, is a classical object, the scale L 
is the only classical one in the problem that should be used in order to make C/h dimensionfree. 
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Coming back to flat space m = (A J a , Ej) = (0, 5?) we want to fix e by requiring that the relative 



fluctuation (|3.8| ) for O = Q{uj) is minimized. This leads to the condition that (notice E = 1) 



e 



3 



[t 2 /L 2 f (3.10) 



Vol(supp(o;)) [e 2 /a 



212 



be minimized where a := L. The fluctuation contribution depends on the volume of the support of 
uo. Since we want to resolve regions with our graph of the linear size bigger or equal than L (think of 
L as the smallest wavelength to be resolved for our applications) we obtain that ( ^.10[ ) is certainly 
dominated by 

t l + le*/L 2 Y (3.11) 

This function has a unique minimum at 



2 1 

P+\ ^ 5 < 2 



?pL 1-a , a = r < ^ < ^ (3.12) 



In |T9| we chose Vol(supp(u;)) > e 3 and different observables, adapted to the graph in question, in 



order to have the lattice degrees of freedom well approximated and led to a ~ 1/6. However, it is 
clear that this choice would lead to boundary effects if the support of u is not adapted to the graph 
in question which would be unnatural. Such boundary effects are avoided by Vol(supp(u;)) > L 3 and 
go at most as the quotient between the volume of a shell of thickness e at the boundary of a region of 
volume L 3 and its volume, that is, as e/L. This drives the lattice scale e closer to the Planck scale. 
Notice that in any case I'p < e < L. 

This concludes the present section. The relations a := L and (|3.12| ) will be our working proposal. 



4 Coherent States Expectation Values 

The purpose of this section is to present the calculation of the expectation values of the various 
terms occurring in the Hamiltonians of section 4 in in the coherent states for QGR discussed in 
the preceeding sections. In the first part we will explain the simplifying assumptions used for the 
computation and introduce the necessary notation. Section |0| is devoted to the computation of the 
expectation values of the volume operator V v and the operator 

Qi(v,r) = ±tr(TMK 1 ,(V v ) r }), (4-1) 



as they are the basic building blocks of the Hamiltonians obtained in |1| . In section [4.3| the results 
are used to give the expectation values of the geometric operators occurring in the Hamiltonians for 
the scalar and the electromagnetic field. 

4.1 Implementation of the Simplifying Assumptions 
The cubic lattice: 

For reasons already explained in our companion paper, the first simplification that we will make 
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concerns the random graphs: In the following we will exclusively work with states based on graphs 
of cubic topology. This simplifies both the notation and the c-number coefficients occurring in the 
Hamiltonians. In a graph of cubic topology, each vertex is six-valent with three edges ingoing and 
three outgoing. We denote the outgoing edges by ej, I = 1,2,3 and choose an ordering, such that 
the tangents of e\, e2, e% form a right handed triple wrt. the given orientation of E. The vertices can 



be labeled by elements v of Z . We write ej 



ei (v — I) 1 where n — I denotes 



the point in Z translated one unit along the negative I axis. In keeping with that convention, we 



associate to e T (v) the dual surface S, 



ei(v-I) 



with its orientation reversed. 



Replacing SU(2) by U(l) 3 : 

We substitute SU(2) by U(l) 3 in our computation because the results of |T(], [TTJ reveal that the 
qualitative features are untouched so nothing conceptually new is learned when doing the much 
harder non-Abelian computation. For the exploratory purposes of this paper it is thus sufficient to 
stick with the Abelian group. Consequently we will replace Q as well as the volume operator itselve 
by appropriate U(l) 3 counterparts. For U(l) 3 each edge is not labelled by a single, non-negative, 
half-integral spin degree of freedom but rather by three integers rij G Z, j = 1, 2, 3 and we have three 
kinds of holonomies h? e . The generators Tj of U(l) 3 are simply i (imaginary unit). The canonical 
commutation relations on L 2 (U(1) 3 , d?nn) are replaced by 

\h j ,h k ] = 
\pj,h k ] =USfhP 
[Pj,Pk] = 

(cf. ( |2.4j) ) with adjointness relations (hpy = {h?)~ l , (pjY = pj. It follows that ( |4.1|) gets replaced by 



Qi(v,r) = -hi 



Finally the expression for the volume operator in our companion paper is replaced by 



Y e t( v ) _ y £ i ^ Y e ^ v) - Y" 2 ^ Y e ^ v) - Y e? ' 

J I r k k I r / / 



with Yf = ih j d/dh j . 

The U(l) 3 coherent states over any graph 7 are given by (see ||10|| ) 



where 



and gi(m) = e p ^ m ''hi(m) EC — {0} = U(l) c . Here m is a point in the gravitational phase space and 



hl(m) = Vexp(i / A 3 ) 
12 



that is, due to the Abelian nature of our simplified gauge group the path system in S e is no longer 
needed. 

As is obvious from the explicit form of the Hamiltonians, our calculation can be done vertex 
by vertex since there is no inter-gravitational interaction between the associated operators. We can 
therefore concentrate on a single vertex for the remainder of this section and drop the label v in what 
follows. 

For the sake of the computation to follow, we introduce the shorthands 

h Jaj = h 3 e „, p Jaj =p/, gjaj = e PJ ^h J(7j 
and similarly the operators Yj a j = Y- J . Let us finally define 

D = ±V, ^(r) = ^«j(r). (4.2) 

Note that q is essentially selfadjoint. 

The huge advantage of U(l) 3 over SU(2) is that the "spin-network functions" 

T {nJa3] {{h Ja] }) = \\h~:r 

Jaj 

are simultaneous eigenf unctions of all the Yj a j with respective eigenvalue nj a j. Even better, the 
operator qj a j ( r ) * s a l so diagonal with eigenvalue 

xr fJ X \ o A li n ^}) ~ Ar (( n ^j + Wojo),(J°j)}) 

A j O a O j y{ n J^i) - 1 ~ t » 

where 



■ n i,+,j ~ n>i,-j, r n 2,+,k - r n 3,+,i - n>3,-,h 



4.2 The Expectation Values of q 

Now we will explicitly calculate the expectation values of powers of the operators q and U . The 
gravitational parts of the matter Hamiltonians constructed in [I] are all sums and products of these 
operators which act only on the edges of a specific vertex, therefore we can restrict consideration to 
a single vertex and consequently to a part 

Jaj 

of the coherent state which just contains the factors corresponding to the edges of a single vertex. 
What we are looking for is the expectation value of an arbitrary polynomial of the q: 

. (^ 3jCTj },nf=i^^ fc (^)^ } ) 
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rwi^ 112 



where (see fllQfl ) 



lWll 2 = \/y / '[l + ^{p)] 1 « = ^ |lf*(p)|<if* = 0(0- ( 4 - 4 ) 



As in ||10|| , in order to extract useful information out of the formula (|4.3| ) it is of outmost importance 



to perform a Poisson transformation on it because we are interested in tiny values of t for which ( |4.3| ) 
converges rather slowly while the transformed series converges rapidly since then t gets replaced by 
l/t. To that end, let us introduce T = ^/t,xj a j = Tnj a j, whereupon 

^' TT Ik/;* II 2 



where 



A J <r j y\ x J^i) - 1 1 



Y({ Xjaj} ) = t»V4 Ai^t MlfW ] [ " W ; X2 '-- fc ] [ ^"^ l (4.6) 



Then Poisson's theorem gives 

TT 11^* II 2 

An observation that reduces the eighteen dimensional integral to a nine dimensional one is that the 
integrand in fl4.7|) only depends on xjj = Xjj = [x — xj-j}/2 and not on Xjj = [xj } +j + xj-j}/2. 
Consider also the analogous quantities p^- = [p ±pj_j]/2, = [n^+jin j_ )3 -]/2 and let pj m = 
Pjj, ^Jm = Pjj- Switching to the coordinates Xjj, noticing that | det (d{xj a j }/d{x~jp x~j-}\ = 2 9 we 
obtain 

(f ) 9 E {wjgj} [/ R9 ^'-W^W] 



X 



n^ii^ji 2 

TV 



/ d V^[-^^^^)/TlTJ A ^ fe . ({lj . }) 
/ R9 fc =i 



(4.8) 



where 



X r ,r X N o ^d^l) ~ A r ({xjj + T^ (Joio ) !(Jj )/2}) 

\ r ({xj J })=t^ 4 (\det({x J ,}) r/2 (4.9) 
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actually no longer depends on cr ' The integral over Xj- in Q4.9Q can be immediately performed by 
using a contour argument with the result 



(4.10) 



n \w ii 2 

1 U,(T,j I l^ffjaj 1 1 

Finally, using (|4.4j) we can further simplify to 

/ - 

^ ; [(1 - JT t ) M , (1 + tf t ) u ] ^ 

r N 

x ^^E^-^^Jife-^ilfflTTy^ (4.11 
Jr 9 fc=1 



where the notation for the denominator means that its value ranges at most in the interval indicated. 
Its precise value will be irrelevant for what follows since its departure from unity is 0(oo). 

Only the nj a ~ = terms matter: 



The remaining integral in ( |4.1 1| ) cannot be computed in closed form so that we must confine ourselves 



to a judicious estimate. We wish to show that the only term in the infinite sum of fl4.11|) which 



contributes corrections to the classical result of finite order in t is the one with nj a j = for all J, cr, j. 
In order to do that, we must demonstrate that all the other terms can be estimated in such a way 
that the series of their estimates converges to an 0(t°°) number. This would be easy if we could 
complete the square in the exponent of the integrand but since for r/2 not being an even positive 
integer the function A r is not analytic in C 9 we cannot immediately use a contour argument in order 
to estimate the remaining integral. In order to proceed and to complete the square anyway we expand 

the product Ylk=i Ki k a k j k ({ x Jj}) monomials of the form Ylk=i A ^ Xj ^ +c Jj^ w ith Cj. = T5j k j k: jj/2 
or Cjj = and estimate the integrals over the latter. We trivially have 

+ c%}) = e/\[det({ Xjj + 4-})] 2 r/ 4 = ^ 4 exp (L ln([det({^ + c^-})] 2 )) (4.12) 

where we must use the branch of the logarithm with \n(z) — ln(|z|) + iip for any complex number 
z = \z\e ltp with ip G [0, 2ir). With this branch understood, in the form ( [4.12 ) the integrand of (|4.11| ) 



becomes univalent on the entire complex manifold C except at the points where det ({ x jj + Cjj}) = 0. 
Now a laborious contour argument can be given tho the extent that we can move the path of 
integration away from the real hyperplane in C 9 without changing the result. Therefore we can 
indeed complete the square in the exponent. 

It remains to estimate Q4.11J) from above. Isolating the term with nj a j = for all J, a, j we have 



■9 

2 „ TV 



< • > " [( i-g,A(i + g,H I n w*« + pv*» 
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E 

{nj<rj}^{a} 
N 



[(l- K t y\ (i + K t ys 

9 fe=i 



;i-A',)is 

AT 



E 



.2 

Jj X 



{"Jaj}^{0} 

x J^^ e iln(|det({Txjj + (pj J -i7rn Jj )})|) g f ln(| det({Txj J +t<5 { j j)i(i7fc ^ ) + (pj J -j7rn Jj )})|)j 



fc=l 



(4.13) 



Let wjj be a matrix of complex numbers and define the norm ||iy|| =Y2jj \ w Jj\ 2 so that i n particular 
||^i + ^2|| < I l w i| I + I |w 2 | | and | to jj | < \\w\\ for all J,j. Now det({u>jj}) is a linear combination of six 
monomials of the form w j^w j 2 j 2 w j 3 j 3 so that | det({wjj})| < 6 \\w\\ 3 . In particular, | det({Txjj + 
(pjj - mnjj)})\ < Q{T\\x\\ + ||p|| + vr \\n\\f and | det({Txjj + t6 {JjU j kjk) /2 + ( Pjj - iimjj)})\ < 
6(T ||a;|| + t + \\p\\ + n \\n\\) 3 . Invoking this result into ( |4.13|) we find 

A 



< - 



< 



t 

A 6 r/2 
" t 



;i-^) 18 



,iV 



E «" 

9 



1 Y*J<Tj n Ja 



(f xe - 2 \\ x \\ e ^r ln(6[r||x||+t+||p||+7r||n||] 3 ) 



:i-^) 18 



E 



-2- V n 2 

g t £-iJrr] '\lcrj x 



{"J CTJ }^{0} 



X 



/" rf 9 xe -2||-H 2 [I +t || x ||2 +t+ || || + 



7r n 



(4.14) 



where [3iVr/2] is the Gauss bracket of a real number (largest integer smaller than or equal to 3Nr/2) 
and in the last step we have used the elementary estimate x < x 2 + 1/4 valid for any real number 
x. The integral in the last line of (|4.14 ) can be evaluated exactly by invoking the binomial theorem. 
Consider the integrals of the form 



h = \l- 

7T 



d m xe- 2 ^ 2 \\x 



i 2k 



(4.15) 



for any positive integer m. Switching to polar coordinates one easily proves the recursion formula 

m + 2(fc-l 



4 



-h-i 



(4.16) 



and since In = 1 we find 



h 
h 



it m even 



2 fc (f)! 

it m odd 



m — 1, 



' to— 1 

. 2 
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+ *)! 



(4.17) 



Using the elementary estimate e(n/e) n < n! < e((n+ l)/e) n+1 we find for < k < n and n > 2 that 

m + 2n /2 m + 2n k . ,m + 2n k 

h < e( 2e ) m/2 ( 4e ) fc = C m ,n( 4e ) if m even 

h < / 2 fry - m 77 TT == ^ n )* if m odd 4.18 

2e (m — lj! m — 1 4{m — 1) 4(m — 1) 



In our case m = 9 and n = + 1. Thus, we can finish the estimate of ( |4.14f ) with 



2 

^9 



'4 6 r / 2 \JV 



(1-i^) 18 ^ 

.1 9 + 2([^] + l) iraWP1+1 

X [ 4 + * 32 + t + M + 71 1' 72 "] ^ 419 ) 

which is obviously of order Oft 00 ). We can give a bound independent of p since in our applications 
||p|| can be bounded by a constant of the order of t a . 

Let us summarize our findings in the form of a theorem. 

Theorem 4.1. Let \\p(v)\\ 2 = PJj( v ) 2 ■ Suppose that there exists a positive constant K such that 
su P«eV(7),me.M =: Ibll < K is uniformly bounded. Then for small t 

-9 

2 „ N 



< • > = [d-x.ra +*-,)'«] /,, ^' 2E "' % n Wfr* wn) + °( t0 °) < 4 - 20 > 

independently of m £ A4,v <E V(-y) . 

Expansion of the remaining integral: 

It remains to compute the power expansion (in T) of the remaining integral in (|4.20|) and to show 
that at each order the remainder is smaller than the given order. We will see that only even powers 
of T contribute so that this expansion is actually an expansion in t. The basic reason is that the 
expansion of the integrand in powers of T is at the same time an expansion in powers of xjj as is 
obvious from the explicit form of the functions X r ({xjj}). These powers of xjj are integrated against 
the Gaussian e -2 "^" which is an even function under the reflection xjj — > —xjj whence the integral 
for odd powers (an odd function under reflection) must vanish. We will not be able to show that the 
integral in (|4.20| ), which certainly converges for any pjj,t (just set ||n|| = in above estimate), can 
be expanded into an infinite series in powers of t, rather our estimates will be only good enough in 
order to show that there is a maximal order no (which becomes infinite as t — > 0) in the sense that 
the remainder at order n is smaller than the given order for all n < uq. We will use rather coarse 
estimates which could possibly be much improved in order to raise the value of hq derived here but 
for all practical purposes the analysis described below will be sufficient since n is anyway a rather 
large positive integer. 
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Consider once more the function X r Ja Ax+p/T): Let us introduce q = pt a which is of order unity 
and s = t 1 / 2 ^". Then 



r/2 \ det(l + q- 1 xs)\ r ^ 2 - | det(l + q~ x xs + q- 1 5 Jj sT/2) 



r/2 



r Jaj (z+p/T) = 2| det(p)r ^ ^ ^ 1 ^ 2 JJ ' " (4.21) 

Now for any matrix A we have det(l + A) = 1 + Tr(A) + |[(Tr(A)) 2 - Tr(A 2 )] + det(A) =: 1 + ^and 
so det( l + A ) 2 = l + 2z' A + (z' A ) 2 =: 1 + z A =■ Va>0. Let y = 1 + Zg-i^ and Z/i = 1 + ^-^xs+^-sT]- 
Then Q4.21|) becomes 

AS^(x + P /T)| = ^M^[^-^] (4.22) 

and we should expand y r ^, y '^ around y = y\ = 1. We now invoke our knowledge that < r < 1 is 
a rational number, so we find positive integers M > L > without common prime factor such that 
r/4 = L/M. Let us define recursively 



fflAv) = v L/M , 

fffiiv) = " !r :\" ' - ■ (4-23) 



f (n+l)/ x . fL/\l(y) /l/A/W 



It follows from this definition that 

f?Uv) = E - if + - !] n+1 - ( 4 - 24 ) 



fc=0 



Lemma 4.1. M^e /iai>e 



/J(l) = (Wi) (4.25) 



where 



(L/M)(L/M- 1)...(L/M -k+1) _ fc+1 L M — L 2M — L (k - l)M - L 



lL , u ,k) = ^-^ ^ '- = (-D' M 2M 3M ... H# 



and the following recursion holds for all n > 1 

V^i-l An) ( \ _ sr^n M) ( i \ v^A/-l f (n-i+1) ,\ 
^(n+1)/ , _ L~ik=l J k/M\y > Z^Z=1 JL/M\ L > Z^k=l Jk/M KM ) ( A OR\ 

J L/M \V) ~ f (0) / x • t 4 ""^ 

2^fc=o h/M\y) 

The proof of the lemma consists in a straightforward Taylor expansion (first part) and an induction 
(second part) and will not be reproduced here. 

The motivation for the derivation of this recursion is that it allows us to estimate \fjjM (y) \ once 
we have an estimate for all the \f^J M (y)\ with < k < M — 1, < Z < n. 
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Lemma 4.2. For all < L < M, n > we have 

\f ( L/M(y)\<(± + y)(m n 

where (3 > 1 is any positive number satisfying (3 > 1 + -Jzi, e.g. (3 = 3. 



(4.27) 



This lemma can be proven by induction, using the results of the previous one. 

Using the expansion ( |4.24| ) and the fact that y is a polynomial in the xjj it is possible evaluate 
the Gaussian integrals over the first n terms the last one of which is obviously at least of order s n . 
We would like to know at which order uq the remaining term in ( |4.24j ) is no longer of order at least 



To that end recall that y = 1 + 2z + z 2 where z = Tr(A) + |[(Tr(A)) 2 - Ti(A 2 )} + det(A) and 
Ajk = s 'Yl l j{<f 1 )jjXjk- We now have the following basic estimates 



r Tr(A)\=s\^2q J jx J A<s\\q x \ 
Jj 



x 



Xjk\ 



Tr (A 2 )\ = s 2 \^2(q l x) jk (q 1 x) kj \ < s 2 \ ^ 10? lx )jk\ \{q l %)kj\ 

jk jk 



Xjk\ 



\ 



^ 2|| — 1 1 1 2 II i|2 

< s \\q || ||x|| 

3 1 1 „ — 1, 



det(A)\ < 6s 3 ||g _1 x|| 3 < 6s 3 ||g 



„3||„-1||3 



X 



where in the first line we have made use of the Cauchy-Schwarz inequality for the inner product 
< x, x' >= x Jj x 'jjy m the second for the inner product < x, x' >= ^2jXjx'j, in the fourth line 
for the inner product < x,x' >= 'Ylj x j x 'j an d finally in the last line we have used the estimate 
derived between equations ( |4.13 ) and (|4. 14 ) . These estimates imply that 



\z\ < s \\q 



-H 



x 



2 II -1|| 2 || ||2 

s \\q \\x\\ 



6|det(g )| \\x 



--: u\ x 



\y-l\< 2u + u 2 =: P(||x| 



and P(||a;||) is a polynomial of sixth order in ||x||. 

We are now ready to estimate the Gaussian integral over the remainder: 



E„ 



7T 

< 



d»xe-^ 2 f^\y)[y-l] 



n+l 



(3M) 



n+l 



d xe 



-2\\x\\ 



[(P(\\x\\)T +2 + 2(P(\\x\ 



\n+ll 



(4.28) 
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Consider an arbitrary polynomial in ||a;|| of the form 

i 

p ( x ) = ^2 a k 

By the multinomial theorem 



\\k 

\x\\ 



n +..+n;=n v u 7 v l/ fe=0 



™fcl ||~||Eife=o fcn * 



Let us consider Gaussian integrals of the form 



/ — m I — m oo / — 171 

^ J^d m xe- 2 ^ 2 \\xf = V m ^l J dre-^r n+m - l =:V m ^ J, 



n+m— 1 



where 1/ m = 27r m / 2 /r(m/2) is the volume of S m . Now 



J n = 3ra/2 ^y for n even, 

4 2 • 



J n = V^" 1 ^^— - -!) for n odd, (4.29) 



and one immediately checks that 



_2^[f] 
4 2lf 



where [.] again denotes the Gauss bracket. Using the above used estimate for the factorial n\ < 
e ( (w+i) ^re+i we ma y f ur th er estimate 



J " " "1 " ~4 



where we used ! ^ — [f ] — f • Finally, if n < n M then 



j, < £^.(=-±1)*. (4.30) 



Combining these results we obtain the final estimate 

-m i — m , i 



— m / — m . i 

2 - [ d m xe-^ 2 p(xT = v m J 2 y w! M <V^ 



no+..+n ; =n v / \ / £,_q 

I 



- ^Vtt 4 2^ (no!)..^,)! 1 !! * J 1 e j 
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V, 



m—1 



— / 

2 ev27r m + /n N 

~1T ^^T' 



i 



k=0 



1+ln 
4e 



,m + Zn. ... 
--:K m A— 



m + In. 



(AM) 



4e V 4e 

since X/l=o — ^ n = n M— m for any configuration of the subject to the constraint n +..+n/ = n. 



In our case we have m = 9, 1 = 6 and thus we can bound the remainder (|4.28 ) from above: 

,9 + 6(n + 2) 9 



E n < K 9t6 (3M) 



n+l 



4e 



rfpf ^9 + 6^+2) )n+2 



(4.32) 



^ 9 + 6(^ + 1) ^^ /l + 6(n + 1) , n+1 



4e 



4 (; 



For small n the error E n is the number s n+1 times a constant of order unity. For large n, however, the 
error becomes comparable to the order of accuracy (in powers of s) that we are interested in itself. 
The value n = n from where onwards it does not make sense any longer to compute corrections can 
be estimated from the condition 



E n +\/E n > 1. 



(4.33) 



Due to the complicated structure of (|4.32|) the precise value of n cannot be computed analytically 
but its order of magnitude can be obtained under the self-consistency assumption that n is quite 
large so that the change of P(y(9 + Q(n + 2))/ (4e)) as we change n by 1 is much smaller than its 
value. A tedious but straightforward estimate shows that under this assumption 



4e( 



n 



MM) \2 



-9 



6 



(4.34) 



where Tq(M) is of order unity. Thus n$ is a very large number if | \q 1 1 1 is of order unity and s is tiny. 
Moreover, 



5P = 2{u + 1)(1 + 2r + 18t 2 )St = 6(u + 1)u8t/t < QP-. 

T 



But under the change 5n = 1 



5t « — 5n = — 

dn 9(9 + 2n) 



whence 



< 



< 1 



(4.35) 



(4.36) 



(4.37) 



V p jn=n _ 3(g + 2no) 

as desired since no is a large number. 

Let us now finally go back to our desired expectation value ([4.20|) which we would like to compute up 
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to some order n < n in s. Let again y = 1 + z q -i xs = 1 + z and yj a j = 1 + z q -i[ xs+ s JjsT / 2 ] = 1 + z j&j 
with z A = (z' A ) 2 + 2z' A , z' A = Tr(A) + ±[(Tr(A)) 2 - Tr(A 2 )} + det(A) for any matrix A and recall our 
convention r/4 = L/M. Thus Q4.22|) becomes up to order n 

X r {x + ,,/T) = 2|det(rt|2L/ " [//«' - (4.38) 



k-l 



2\det(q)\ 2L / M t 6L / Ma „. ,^Ak) \ - 

k=l 1=0 

+ [ff/lP(y)(y - i) n+1 - /;.";; - i) rt+1 ]} 



k-l-h 



In order to compute ( |4.20|) up to order n with respect to s we have to plug the expansions (|4.38|) into 



formula ( |4.20| ) and to collect all the contributions up to order s n . The integral over the remainder 



is then still smaller as long as n < tiq as shown above. In the present work we are interested only in 
the leading order (classical limit) and next to leading order (first quantum correction) as well as in 
an estimate of the error at the next to leading order. 

A laborious but straightforward power counting reveals that 

X r Jaj = ^{l + sx + {sx) 2 + 0{sT)) (4.39) 

where the notation just means that Aj CTJ - is a polynomial in xjj of order two where the monoms of 
order 0, 1, 2 come with a power of s of the order indicated or higher. We thus see that if we wish to 
keep only terms up to order (sT/t) N and {sT/t) N s 2 in Yik=i Kj k a k j k (. x +P/T) it will be sufficient to 
do the following: For the term of order (sT/t) N keep only the terms proportional to x° in each of the 
factors of the form fl4.39|) and for term of order (sT/t) N s 2 keep either 1. only the terms proportional 



to x 2 in one of the factors of the form ( 4.39 ) and only the terms of order x° in the others or 2. only 
the terms proportional to x 1 in two of the factors of the form ( |4.39| ) and only the terms of order x° 
in the others. Clearly terms of order (sT/t) N s do not survive since they are linear in x and integrate 
to zero against the Gaussian. 

In estimating the error that we make notice that there are two errors, one coming from neglecting 
all higher orders in ( |4.39|) and one from the remainder in the expansion (|4.38|) . As for the first error, 



notice that all Gaussian integrals are of order unity so that the indicated powers of t correctly display 
the error (compared to (sT/t) N s 2 ) as of higher order in s. As for the second error we can use the 
expansion (|4.38|) up to some order n! > 2 until s n+l sTs 2 in view of the estimate d4.32|) . The 



minimal value of n' depends on the value of a. For instance, if a = 1/6 as indicated by [IS] then 
s = t 1 / 3 so that s n '~ 2 = t( n '~ 2 )/ 3 « T = t 1 / 2 means n' > 2 + 3/2 so the minimal value would be 
n' = 4 in this case. This value is well below n ^> 1 so that the error is indeed of higher order in s 
as compared to (sT/t) N s 2 . 

With these things said we can now actually compute the first contributing correction to the 
classical limit. We will not bother with the higher order corrections since we just showed that they 
can be bounded by terms of sub-leading order as compared to (sT/t) N s 2 . In particular, we will 
replace the 0(t°°) corrections by zero in ( 4.20|) . We then have 



— 9 , N 
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N 

z=l fc/« 
+ [ £ ^(^^^lA^^J.+p/T),^] J] A^(, + p/r) |l0 ] (4.40) 

+ 0(t (Ar[3r/2 ~ 1]a sT) (4.41) 

where the restrictions mean the ones to the appropriate powers of x as derived above. It remains to 
explicitly compute these restrictions and to do the Gaussian integrals. According to what we have 
said above we write 

AS CTi (* + V/T) = 0(t^-^sT) + 2| det(g)r/ 2 t^ 2 -^{[/ r ( !)(l)(^|^) |:c o] 



Jaj\ ' irl J *s yv o± j -r i| <acL V L/;| v l[J r /4V- L A TjT 

+ [/$(i)(^f^M + / r ( / 2 i(i)(^| M )|,o((i/ - 1)1.1 + (y^- - i)|.01 

+ [/$(i)( 2 4r^M + /^ax^^MO/ - + - 1)1--) 

+ Oi)^)^^ - i)i*0 2 + (fe - i)i«0 a + (y- i)\Avj« - i)i,0]} 



(4.42) 

And furthermore 

y - 1 = 2sq^ m x Mm + s 2 {2q^ m q^ n - ql^q^JxMmXNn + 0(s 3 ) 
=: sC Mm x Mm + s 2 C Mm ' Nn x Mm x Nn + 0(s 3 ) 
j/jaj - 1 = 2sTr(q~ 1 x) + s 2 [2Ti{q- 1 x) 2 - Ti(q~ l xq~ l x)} + 0{sT) 
=: sC Mm x Mm + s 2 C Mm > Nn x Mm x Nn + O(sT) 

2 

Vj03 sT V = Qjj + S ( 2 ?J"/?Mm - Vjm<htj) x Mm + y {^{q- l )e jrnn ejMN + Qjj {QlrfmQNn - VmuQnU 

iNn ~ Qj n lNj)] X Mm x Nn 
= '■ Cj a j + sCy™XMm + s2 C J( ™' n XMmXNn- (4.43) 

We can therefore simplify (|4.42 ) to 

X r Jaj (x+p/T) = 0(t^M a sT) + 2| det{q)Y^I 2 -^{{fll{l)C, J(TJ ] 
+ s\f%Q)C%? + 2f^(l)C J(7J C Mm ]x Mm 

+ s 2 [f^l(l)CZT Nn + 2f^{l)C, Ja] C M ^ Nn + 3f^(l)C J(TJ C Mm C Nn }x Mm x Nn } 

=: 0(f r ' 2 ^ a sT) + 2| det(q)\ r / 2 t^ 2 ^ a {D JaJ (r) + sDf a J(r)x Mm 

+ s 2 D^ Nn (r)x Mm x Nn }. (4.44) 

Putting everything together now yields the following theorem. 
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Theorem 4.2. For the classical limit and lowest order quantum corrections of expectation values of 
monomials of the operators qj a j (r) for topologically cubic graphs we have 

N 2 N 

k=l M,m 1=1 k^l 

+ E ^WCMll^W]} (4-45) 

\<i<l<N k^l,i 

where the constants are given by 

r Mm,Nn _ r, -1 -1 _ -1 -1 
Cjcrj = Qjj ) 



^Jaj 


= i[det(g 


D J*j 




= fr/i(^)Cjaj, 


r\Mm 
U Jaj 






Mm,Nn 
Jaj 


J") 


/■(l) /i \ n Mm,Nn 



ml 



^QMrni Ijj iNn 1 Jn In] ) 1 > 



f(2) /-r \ /-» riMm 

D2?™(r) = f%{l)C%?" n + 2tfl(l)C Jaj C Mm > Nn + 3f^(l)C Jaj C Mm C 
and the fj!j\(X) = (r/4, fc) are simply the binomial coefficients. 

The first correction is small as long as a < 1/2. T7ie error as compared to the first quantum correction 
of order 0(t^ N ^ 3r ^ 2 ~^ a s 2 ) is a constant of order unity times f( N l 3r / 2 - 1 ] a s T and thus small as long as 
< a. 

So far we did not look at the classical limit and the first quantum corrections of (powers of) the 
volume operator itself but it is clear that it can be analyzed by similar methods, in fact, the analysis 
is even much simpler because we just need to expand X r (x +p/T) in powers of s without dividing 
by t and thus will have to do an expansion in terms of y — 1 of one order less than for X r Ja -{x-\-p/T). 
Clearly the classical order will be that of |det(p)| r / 2 = | det(g) | r / 2 f 3 ™/ 2 = 0{t 3ra / 2 ) and the first 
quantum correction will be of order 0(t 3ra ^ 2 s 2 ). We thus have, in expanding up to second order in 
y — 1, where y = det(l + sq~ 1 x) 2 as before 

Y(x + p/T) = |det(g)r/¥~/ 2 |l + S / r ( / 1 i(l)C' Mm a ;Mm (4.46) 
+ s 2 [tf\(l)C Mm ' Nn + ${(l)C Mm C Nn ]x Mm x Nn ) + 0(t 3ra ' 2 s 3 ). (4.47) 



Thus we obtain an analogue of theorem 4.2 above 
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Theorem 4.3. For the classical limit and lowest order quantum corrections of expectation values of 
powers of the volume operators U v for topologically cubic graphs we have 



(ib U ib ) 

% U2 = \det(q)\ r/2 t 3ra / 2 {l + j Y,lfr%( 1 )C Mm ' Nn + fl%{l)C Mm C Nn }. (4.48) 



,2 



} ' ' M,m 



The first correction is small as long as a < 1/2. The error as compared to the first quantum correction 
of order 0(t ( - N ^ r ^ 2 ~ 1 ^ a s 2 ) is a constant of order unity times t^ 3r ( 2 ~^ a s 3 and thus small as long as 
< a. 



4.3 Application to the Hamiltonians 

So far our considerations were completely general and model independent and we see that our coherent 
states indeed predict small quantum predictions as long as < a < 1/2 and lp/L <C 1 with 
controllable error. However, now we will specialize to the case of scalar, electromagnetic and fermionic 
fields coupled to gravity and compute the expectation values of the relevant gravitational operators. 

We recall from our companion paper |l[ that on a cubic graph, the effective Hamiltonians for the 
scalar and the electromagnetic field are 



cff 
KG 



2Q 



|-E (<^>)>^- E ( 9 t^U v ) (d+,\nU v )-K 2 (V v )(lnU v 

v \ Ial'a' 



eff 
EM 



2Q 



- <^f>)>E hh^v)] [E l < H ?i W2 m h 



(4.49) 



(4.50) 



where ( ■ ) denotes the expectation value in a semiclassical state for the gravitational sector, and the 
geometric operators are given by 



FkM = TTTT 



1 , 



1 , 



i j2 ^^^^(t,,;)^^;)^^;) 



der 



mag \ u I 



--Y 



fl,(T2,(r3 

UK 



0"2,0"3 



1 1 

i i 



64 i% 



E#>^)#'>^). 

3 



The matter fields are represented as 



tt v = ihQ K GY v , E Ia (v) = ih~Q EM Y Ia , B Ia . aulT3 ( 
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where the Y are invariant derivatives on U(l), U v is a U(l) point-holonomy and Hp a U(l) holonomy 
around a minimal loop. 

The Hamiltonian for a fermionic field is given by 



k ^ = ~W £ [Ob(v')§\(v) - B' B (v')9'l(v)) x 



x {\tiike IJK E ^/oi( u ^)^( u .5)[ rfc (^( u )^,/(4 s W)-^,«)]Afl} 

0"1,0"2,0"3 

-{^^ E [([^V)r\ /(e 4 K)) - («', ^)}} 

-ihK £ W«y[^(«')^(«)-^fl(«')^(«)]- ( 4 - 51 ) 

«,«'6V(7) 

We strongly recommend to take a look at ]I] where the above Hamiltonians are derived and all the 
ingredients are defined and discussed in detail! 

We now proceed to compute the expectation values of the geometric operators in a coherent state. 



To this end, we will use the formulae given in theorems [4.2| and [4.3| with the appropriate values of 



r, N, Jk, o~k,jk inserted, and and perform the additional computations necessary. 
4.3.1 The Kinetic Term 



For Fki n we have to use theorem L2 with N = Q. Employing the relation (|4.2p between q and Q we 
find 

1 1 /2ta 3rx6 



<*idn> = jvjj^ ) (2| detCg)! 1 /^^ 1 ^) 6 e J ^e jd; ■ ^ 



6 2 6 



x {[\[Dj k „ kh {vn + j EE<S Mm (V2) n^ fc (i/2)) 

fe=l Af,m Z=l 

+ E ^(V2)^5,(V2) II ^; fe (V2)]}. (4.52) 

l<i</<6 kjfeZ.i 



For r = 1/2 we have 

= /«(!) = i, 02 = /$(!) = -II = 03 = /«(!) = ^ = _|_ (4.53) 

and consequently 

E <? Mm (l/2) = [O! + 3a 3 ]g7/ Tr^ 2 ) - ^ (4.54) 
E ^• 1 (l/2)^£ J ' 2 (l/2) = 4[ai + a 2 }%\ iq j\ 2 Tr(g~ 2 ) 

Mm 

- 2a 1 [a 1 + a^qj^qjX + q^L,:) : ) + a\qf xJ ^ h . (4.55) 
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Now we have to deal with the contractions in ( 4.52 ). It is easy to see that 

6 



..J1J2J3 J4J5J6 



36 



e e 3ihh e e 34jde 



IK', 

fc=i 



J4J " 6 fe=i ' det(g) 2 ' 

12Tr(g" 2 ) 



e JlJ2j3 W 3 e J4J5J6 W 6 ^^) II = if i,i 6 {1,2,3} or J, i G {4,5,6}, 

4Tr(g- 2 ^ 



J1J2J3, ,^4J5^6, „— 2 -2 

fc t jlj2J3 t t 3i3S36 L LJ i .J l ( lj i j l 



n si 



det(g) s 



otherwise. 



(4.56) 



Using the above together with ( 4.54]) in Q4.52|) yields 



(Fkin(«)) 



4 6 6 (2|det(g)| 1 / 4 t[ 3 / 4 - 1 ] 



det(g) : 



{ 36[ a ?] + ^T r(9 -)K(36[a 1 + 3 a3 ]-12| 



+ a 4 (15(4[ Gl + a 2 ] 2 (36) - 2a x [ai + a 2 ](12 + 12)) + 9a 2 ]} 

^^={1 + l - Tr(jT 2 )[(5 + 24a 3 ) + 15(4[1 + 8a 2 ] 2 - |[1 + 8a 2 ]) + 1]} 

a 9 t 6 1 r 1707^, 

{1+t— Tr(|T 2 )}. 



Let us finally transform back to the dimensionfull quantity P = a 2 p. We get 



y/det P{ 



^ 1707 2 

1 H it F (v) 

t 512 v ; 



(4.57) 



4.3.2 The Derivative Term 



The derivative term F^er requires N = 4. From theorem 4.2 we find 



(Fl: J 'n = ~ ( i 



aa ' 1 (i) 4 ^^!^^! 378 ^ 11 ") 4 E^S^^S^ E ( 4 - 58 ) 



4 4 V3 

x {[II^, fc (3/4)] + J EECf m ( 3 / 4 )II^(3/4)) 

fc=l M,m i=l fc^Z 

+ E ^(3/4)^,(3/4) J! ^,(3/4)]} 

l<i</<4 fe^/,i 



(4.59) 



For r = 3/4 we have 

^„(i )m _3 ^,(2) m 3 29 3-29 ^ f (3) m _ 3 • 2945 _ 3 2 - 5 • 29 



)13 



(4.60) 
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Furthermore the reader may verify that 

4 . 2 

E,JJiJ 2 , ,J'J 3 J 4 , TT „ -i i _ ^IjJ' 
i fc=i w 

Y^ f jJiJ 2f ^'j 3 j 4 ^ _-a ^ TT„-i _ 2 [gjj' Tr (g~ 2 ) njjjj 

Z_> fc iiii2 fc t jj3j4 ( dj l j l ) liyjfci* det(g) 2 ' 

E eJJlJ2e ^ 2 e J ' j3j4 ^ 4 g7i^ II = o if M e {1,2} or /• i c {3,4}, 

i MM 
£ £ « \,/ 4 S.*^%l II l& = & ' T ^I +<W "ttewto. (4.61) 

Thus we can finish with a tedious but straightforward computation: 

(Fl° J '°') = ~ (0 4 ^(2|det(g)| 3 / 8 t[ 9 / 8 -^) 4 x (4.62) 
x fb 4 4g ' J/ 1 I ^l^dn I ^ 1 4q2jJ ' Tr(n- 2 ) «l 2[g|^ Tr(g~ 2 ) - ^] 

x \ K d^ ] + T [4Giaai + 3a3j d^F Tr(g )_ T di%p } 



+ af(4[ai + a 2 ] 2 Q ^ Tr(g" 2 ) - 2a x [a x + a 2 ]12 



j j' 



det(g) 2 det(g) 2 
det(g) T 



2 g^Tr(g- 2 )+^ 



era' t 4 a 9 1 

X 



4 $ ^detM 

x {pjj' + {\Pjj> Tr(p- 2 )[4(1 + 16a 3 ) - ~ + jj(3 + 16a 2 ) 2 - 4(3 + 16a 2 ) + 1] 

+ <^ + 4(3+16a 2 ) + l]]} 
<j(t' t 4 a 9 1 



- 4 4 vraK^&'^^lH}- <« 3 > 

Again as a last step we transform to the dimensionfull quantity P = a 2 p: 

(Fj:/°) = — , 1 {Ph + ^ \—ph> MP' 2 ) + ~8jA I- 

4 ^|det(P)| I * L 128 JJ V ; 32 JJ J J 

4.3.3 The Mass Term: 

We now consider the mass term. Its basic building block is the volume operator itself, so we can 
apply theorem |4.3| with r = 1 . In the by now familiar way we find 



4 

Af.m 
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a 3 |rfet(g)|V^/ 2{1 + ^ Trr2)[ l_ 4 | ]} 



a 3 |cfet(p)| 1/2 {l 



y/detP(v) 



1 + 



|Tr(^) } 
i\ 1 



TrP~ 2 (i; 



(4.64) 



4.3.4 The Maxwell Hamiltonian: 



The operators F e i and F mag differ by their c-number coefficients, but the gravitational operator at the 
heart of both is the same, corresponding to iV = 2 and r = 1/2. In both cases we have to compute 
(^•(1/2)^,(1/2)). 

Let us use the definitions of ai, a 2 , a 3 given in ( f4.53|) and equations (|4.54|) , (|4.55|) . We find 



< ^,(1/2)^,(1/2) >= <^,(2| det(g)| 1 ^[3/4-i] a)2 x 

n ^wi/2)] + j e [E<^ Mm (V2) n^ fe (V2) 

fc=l M,m i=l fc^i 

+ E ^.(1/2)^,(1/2) J! ^*(V2)] 

l<i<K2 fc^,j 

(2 0l | det(g)|^[ 3 / 4 -^) 2 {g 7i 2 72 + £[2([1 + 3^]^ Tr(g-) - ^ 



4[l + ^ 2 + ^ 2 Tr(g- 2 )]} 



(I det(g)| 1 ^[3/4-i] 74 )2 {g7i2j2 + ^[ ? -2 j2 Tr (g- 2 )(7 + 3 
__I^Ml {p-2 Ja + tf^gj^ Tr(p- 2 ) - ^Pj^J}. 



35 7 3 2 • 5 2 



2 7 2 



)- ?^(5 -t)]} 



26 7 ^JlJ2\ 



(4.65) 



We can now employ this result to give the explicit expressions for (F e i) and (F e i). Upon using the 
above expectation value, we find that 



. , era' r r, — ^, ; ^ 9 /763 

w ) = — [v^pm^ + f (512 



Pr 2 Tr p- 2 - —Pr,f 



(plal'a'K 
\ mag / 



64 



V / ^O^P- + ^f^P-TrP- 2 



t V512 



13 

" 16 
13 

16 



p-4 



4.3.5 The Fermionic Hamiltonian 



Due to explicit dependence of (|4.51|) on h' e the expectation values computed so far are not quite 
sufficient in order to compute the full expectation value of the fermionic Hamiltonian. Fortunately, 
the Abelian nature of U(l) 3 allows for a simple transcription of theorem 41 to this slightly more 



complicated situation. Notice that at this point coherent states are, for the first time, essential, 
because weave states, being spin-network states would result in zero expectation values. 
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Theorem 4.4. For the classical limit and lowest order quantum corrections of expectation values of 
monomials of the operators qj a j{ r ) times a holonomy operator for topologically cubic graphs we have 



{gjyjV Jomm Iifc=i HJ k a k 3k\ k '^{g,j a j} _ e -t/Ayi det(g)| r//2 t' 3r//2_1 ' a ) Ar x 



N 2 N 

x {[\[Dj k „ kjk {r)] + ^ ^E^WII^W) 

k=l M,m 1=1 k^l 

+ E D tZn (r)DZn (0 II D ^ Ml W^W/< 

l<i<KJV fc^Z.i 



- ^{ gJgj }>rifc=lgJfc<Tfcjfc( r fc) fe J og ojo^{ gjg; ,-} > _ -t/4iu (r,\ x|r/2.[3r/2-l]miV v 

AT 2 N 

fe=l M,m 2=1 fc^Z 

+ E %(K,wlI^( r )]W«M< ( 4 - 66 ) 

where the constants Dj a j(r), Dj^(r), Dj I J! j L,Nn (r) are defined in theorem^T^ while f^ 4 (l) = (r/4,k) 
is simply the binomial coefficients. The first correction is small as long as a < 1/2. The error as 
compared to the first quantum correction of order 0(t^ Ar ' 3r//2_1 ' a s 2 ) is a constant of order unity times 

£(iV[3r/2-l]a s J> an( J sma H as \ on g as Q < a . 

Of course, in computing the quantum correction in terms of p or q rather than p' = p±iMj 5j j t/4 
or q' = p't~ a up to order t or s 2 respectively one is supposed to insert this substitution into Q4.66Q 
and to drop all higher order terms. 
Proof of Theorem |4.4| : 
We begin with the operator identity |§ 

and exploit that our coherent states are eigenstates for #j 0(7W0 with eigenvalue gj oao j - Moreover, our 
coherent states are expanded in terms of momentum operator eigenf unctions on which the operators 
g-MPJo^ojo and q.j k o k j k (jk) are simultaneously diagonal. It follows that 

N N 



k=l k=l 
N N 



k=l k=l 

It is therefore sufficient to consider the expectation values 

N 



k=i 
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where in the second step we have again performed a Poisson transformation with periodicity param- 
eter T = \/t (see the companion paper for more details). 

As in U we introduce the coordinates Xjj = ± xj-J)/2 and similar for Pj^rij- with 

Xjj := Xjj,pjj := pjj,rijj := nj-. Then one can split the eighteen dimensional integral into two nine 
dimensional ones with the result 



< . > 



t 9 M Ml 

llr {9J<Tj} M n Jrjj 



J2if t^x + e 2E « [ ~ ( *^ )2+to J>Ji~ im Ji + ^ iJ <>* /4)/:zl ] x 



- N 

x [ / ^xe^^t-^+^^ib^-^^+^o^vo.o/^/T] -Q (xj./t)). (4.70) 

k=i 



Then, using the norm of our coherent states and dropping the 0(t°°) terms which come from the 
ones with Ylijaj n2 Jaj > we find 

I — 9 

< . > = J! e f E Jj [(P + +^Joio*/4)^ + (p+^otf Jojo t/4)5 J --(p + )5 j -^] x 

x [ f *h J ! AJfc (x + P+yg ^ /4 )] 

fc=l 

= Jl e^o^oe^l I d 9 xe~ 2 ^ *h f[ X^Jx + P + ^^ )]. (4.71) 



fe=l 



Combining (|4.68|) , ( |4.71 ) we see that compared to (|4.45|) and the prefactor of e t ^hj oao j the remain- 
ing integral in (|4.71|) is the one for the expectation value of the operator monomial Ylk=i QJk<Tkjk( r k) 
just that we have to evaluate it at p + vcr tdj j /4: instead of at p. 
□ 



We are now ready to apply theorems |4.1| and |4.4| to the case at hand. For r = 1/2 we have 

ai : = /«(!) = I a 2 := = --- = - — , a 3 : = /< 8 >(1) = — - = — (4.72) 

1 Jl / 8K ' 8 1/8 8 16 128 1/8V ' 128 24 1024 v ; 



and 



D Ja] {l/2)=a l q- J j 

Df™(l/2) = a 1 (2qj}q^ m - qj^qj*) + 2a 2 qjjq^ m 
= 2[ai + a 2 \q-)qll m - atfj^qj^ 
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= 17 det(g~ 1 )e imra e JM jv + [-^- + 2a 2 + 3a 3 ]g7/gj^g^ 



^ , ~JU0IV~U 1V1 IV 1 L 2 ' ^ ' ^J^JJ 

y + 2a 2 \q Jj q Mn q nm 



- + ^qj-QMnQMrn ~ ^ <i Jn<i Mm<i Nj ■ ( 4 - 73 ) 



It follows that 



Er,Mm,Mrn/i / \ r , Q i -1 -1 -1 a l -1 -1 -1 
^ (1/2) = [<H + 303]?^- q M m<lM m - -j q J^ q MmQMj 



M,m 



[ai + 3a 3 ]gj/tr(g 2 ) - ^g^ 3 



2 

E L, & 1 (V2)^ 2 (l/2) = (2K + a 2 ]q^ n q M l m - a 1 g7 i 1 m g^ l )(2[a 1 + afe 1 ^ - a^q^ 



Mm 
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,.-1 „-2\ 

1 „-3 1 „-l „-3 \ 1 „2 -2 -2 



4[ai + a2] 2 g Jl 1 il g j2 1 i2 tr (9 



- 2ai[ai + a 2 ](g7i 1 ?74 + ^i^J + (4-74) 
The relevant quantity for the Dirac Hamiltonian is 

»,»'£V(7) 



X {< [8v>,e°,{v,l)(cTj[h + hjak ^ T k ]) AB - 5„',e5(„,0)^As]?Mm(^, l/2)<?jVnK 1/2) > 

- < qMmW, l/2)q Nn (v', 1/2) [5^(^,1) ([1 2 + hjak ^ V ^ -T^a^AB - 5 Vje * {v , fl )8 A B\ >} 

+ 2k 5 AB 5 vy [e' B (v')(e A (vy - e B (v')(e' A (v)y}} (4.75) 

or more explicitly 

^6 - y E [-2fcoW„ lW ' feKXM^ - fe(«0 A(«)) f ] 

D,D'sy(7) 

= -f E ^ MJV Wy[^K)(^W) t -^(^)(^W) t ]E x 

u,t) / ev(7) "' ct 

X {[(5„',e5( V ,l) ( <7 j[ 1 2 + -^-^D^B - Sv',e°j(v,0)8AB) < qMm(v, l/2)q Nn (v, 1/2) > 

+ {bvi,e° J {v^){v j Y j Tk\)AB < h Ja k {v)q Mm (v, l/2)q Nn (v, 1 /2) >] 

-[(<^,e3(u',l)([l2 + -y-7" fe ](7j)AB - S^^^Sab) < ?Mm(l>', l/2)q Nn (v', 1/2) > 

+8v,e°j(v>,i)(-^:T k (Tj)AB < qumip' , l/2)q Nn (v', 1 /2)h J(jk (v')~ 1 >]}. (4.76) 
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Let us write q = p(v)t a , q' = p{v')t a , q x = p\{v)t a ,q[ = pi{v')t a where P\{y) = p{v) + 
at5jk/4:, p\{v') = p{y') — at5j k /^- Then, using theorem fO| , we may write ( |4.76| ) in the reduced form 

v,v' £V(-f) 

= -y £ ^e rmny [4(^)(^M) t -^K)(^(^)) t ]E x 

»,«'ev(7) o" 

X {[(^v',e^{v,l)(o-j[l 2 + —^-T h \) AB - 5 v/ ^ vfi) 5 AB ) < g A fm(l/2)gAr„(l/2) > q 
+e~ t/i (5 v '^ j(Vyl) (aj hja ^ Tk ]) AB < q Mm {l/2)q Nn {l/2) > qi ] 

— [(<5i>,e}(V,l)([l2 + - ^ - - Tk]&j)AB ~ ^v,e"j{v' fi)^ Ab) < 9Mm(l/2)?Afn(l/2) >q' 



+e- f/4 ^, e}( ,M) ( 1 r k a 3 ) AB < q Mm {l/2)q Nn {l/2) >,,]}. (4.77) 



It remains to apply theorem [4.4| . We have explicitly for arbitrary invertible q 

e JJ^ 6 mh < g Jm (l/2)g j2J1 (l/2) >= e ^V^ 2 < g JlfflJ1 (l/2)g Jra (l/2) > 



= e JJlJ2 e^ 2 (2|det(g)| 1/4 t [3/4 - llQ ) 2 x 

x {[n^.,(i/2)] + ^ EEcriv^ n^^(V2)) 

fc=l M,m 1=1 k^l 

+ £ ^(1/2)^,(1/2) II ^u fe (l/2)]} 

l<«</<2 jfc^,i 

= e JJlJ2 e^ 2 (2|det(g)| 1/4 t [3/4 - 1]Q ) 2 x 
2 

x {aiffJi^Ji + ^[^(^(K + Sa^tr^ 2 ) - + g ^([ 0l + Sasfe W<T 2 ) - ^J) 

+ 4[ ai + a 2 ] 2 g7 i 1 Ji g7 2 1 J2 tr(g- 2 ) - 2a 1 [a 1 + ^(g^g^ + g^g^J + a??* 2 ^]}- (4-78) 

We have 

J^jMqj^qj^ = 2det(g~ 1 )g Jj 

e JJlJ2 ^ li2 g7i^7i = det(g~ 1 )(g J ,tr(g- 2 -g7/) 

e JJlJ V^ 2 g7j 2 = 0. (4.79) 

Thus we can finish ( [4.78|) with 



e JJ ^e"^ < q Jdl (l/2)q j2jl (l/2) > 
(2|det(g)| 1/4 t [3/4 - 1]a ) 2 det(g- 1 ) x 



,2 



x {2a 2 gjj + -[2a x {2[a x + 3a 3 } qjj tr(q- 2 ) - -[ gjj tr(g- 2 ) - qj}]) 
+ 8[ai + a 2 ] 2 gjjtr(g~ 2 ) - 4a x [ai + a 2 ][gjjtr(g~ 2 ) - qjj}}} 
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4r a / 2 /^/d^){2a 2 1 qj j 

s 2 

+—[(4a 1 [a 1 + 3a 3 ] + 8[<n + a 2 \ 2 )q Jj tr(q 2 ) - {a\ + 4a x [ai + a 2 ])(g Ji tr(g 2 ) - g, 7 /)]} 
+ 24a 3 ] + 8[1 + 8a 2 ] 2 ) gjj tr(g- 2 ) - (1 + 4[1 + %a 2 ]){qj 3 ti{q~ 2 ) - g^ 1 )]} 

/--a/2 2 o2 . e: o4 10 

-.{2 qjj + -[(4 + -g-] + -,)g Jj tr(g- 2 ) - -fetr^ 2 ) - Ql})}) 



I6y/M(g) J 4 LV 2 5 J 2 5 ^ JJ ^ ' 4 
={PJi + 7^[75^tr(p- 2 ) + 52pj})]} 



r Q / 2 , s 



2 



Ujj + ^[75 gji tr(g- 2 ) + 52gT 1 )]}. (4.80) 



8v/d^) 2 7 ' 

Notice that the classical limit is precisely the correct one while the relative first quantum correction 
is given by approximately 1.0s 2 Sjj for flat initial data. 

Now we should compute the additional corrections arising when expanding 



-t/4 JJiJ 2f jjlj2 



< 



9Jiii(l/ 2 )gj 2 j 1 (l/2) > q -* q+ tl-°>va 5j o:io /4 (4-81) 



at q up to order s 2 . However, it is clear that the the additional correction in e~*/ 4 — 1 = s 2 0(t 2a ) 
and the one from 5q = s 2 0(t a ) are both of higher order in s so that we can drop the factors of e~*/ 4 
and the substitutions q — > q%, q' — > q' x in ( |4.77|) which therefore can be written, up to order s 2 as 

v,v'£V(-f) 

X {[{5 v >,e< T J {v,l)(&j[)-2 H r T k ]) AB — SvI^v^Sab) < ?Mm(l/2)<ZiVn(l/2) > g=g(«)] 

— [(^,e5(«',l)([l2 H ^ T k }aj)AB — Sv^jW^Ab) < ?Mm(l/2)gAT n (l/2) >q=q(y')]} 

= -f E fc(«o(^(«)) t -fe(«o(to) t ]E x ^ 4 - 82 ) 

u,-u'eV(7) o" 

x {[(^',e5(«,l)( cr i[l2 H ~~77T" Tfe])AB ~ ^u'.e^u.O^AB) 



x(^- + ^PjMp' 2 ) + 52pj/)])(u)] 
r/x /n , hjvkiv')' 1 - 1 1 

-[( d «,e5(«',l)U 1 2 H 2^ T k\Vj)AB - O v ^ j<yVlfi) d A B)^-J^= X 

x(p/i + ^pjMp' 2 ) + ^p'j])W)}}- 
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5 Towards Dispersion Relations 



In the present section, we will bring together some of the results of the companion paper and the 
previous section: We will compute corrections to the standard dispersion relations for the scalar 
and the electromagnetic field resulting from its coupling to QGR. The necessary calculations are 



performed in section |5.1| for the scalar and in |5.2| for the electromagnetic field. Similar computations 
can be performed for the fermions but they give no new insights so that we leave this to the interested 
reader. We have set up the problem in such a way that the calculations for an arbitrary background 



metric but for a start we confine ourselves to the flat one. In section 5.3 we will discuss the results and 



compare them to those obtained in |25|, |26[. In our companion paper we have given some conceptual 
discussion of the issues involved in obtaining dispersion relation from QGR, so we will mainly focus 
on the concrete calculations. 

In |lj we have obtained Hamiltonian operators for the matter fields of the form 

**f = \ E Hv)P U \v,v')fr{v') + Uv)Q U \v,v>)Uv'), 

v,v',l,l' 

where the coefficients P, Q are the expectation values of specific operators on the gravitational Hilbert 
space. We have computed these expectation values in the preceeding section. 

Note that these Hamiltonians are normal ordered with respect to the annihilation and creation 
operators defined in [[I]]. Thus, the expectation value of these Hamiltonians in a coherent state 
peaked at a specific classical field configuration will yield precisely its classical value. Therefore, in 
discussing the dispersion relations, we will assume the matter quantum fields to be in a coherent 
state and can effectively work with the classical fields p, q. A similar argument can be given for the 
fully quantized Hamiltonians of flU, only that one has to consider a coherent state for the combined 
system of quantum matter and quantum gravity displayed in section 4 of |]] as well. 

Summing up, in the following we will investigate Hamiltonians of the form 

= \ E Vi{v)P ll \v,v') Vv {v') + qi {v)Q u \v,v')q v {v'). (5.1) 

v,v',l,V 

The coefficients P, Q can in principle be taken to be expectation values in a coherent state for the 
gravitational field peaked at an arbitrary point of the classical phase space. However, since we are 
interested in dispersion relations, a notion that by definition describes the propagation of fields in 
flat space, we will restrict considerations to the case of GCS approximating flat Euclidean space 
(denoted by ^fjat in the following. Also, when considering application to situations such as the 7- 
ray burst effect, the curvature radius is always huge compared to Planck length and does therefore 
not lead to any new quantum effects but just to classical redshifts which can easily be accounted for. 

Let us choose the canonical Euclidean coordinate system as global coordinates on E. In the 
U(l) 3 -setting, we can model the flat space situation by choosing the classical values 

Ai(x) = 0, E a T {x) = 51 for all x G S 

with respect to our global coordinates. Therefore all holonomies are trivial and for the fluxes we find 

P<(f) = 4 / dni - 
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We will also use the dimensionfull quantity -Pf(f) = a 2 p\{v ). 

Let us come back to the discussion of ( |5.1| ): Since the coefficients in these Hamiltonians vary 
from vertex to vertex, the equations of motion induced by ( |5.1|) are still highly complicated and an 
exact analytical treatment is beyond the scope of the present paper. Moreover, the solutions to the 
equations of motion will not have the character of plane waves, so the notion of a dispersion relation 
is ill defined anyway. 

In [p]] we argued that in the limit of low energies or, equivalently, large wavelength the field 
propagation induced by (|5.1| ) can be described by a dispersion relation: The graph 7, the GCS is 
based on, breaks Euclidean invariance. However on large scales this invariance is approximately 
restored. 

As we can not easily Q5.1|) compute the solutions to the equations of motion of ( |5.1|) and show that 
they reduce to approximate plane waves with a specific dispersion relation in the low energy limit, 
the question is how one can nevertheless obtain the dispersion relation governing the propagation for 
low energies. 

In |I| we have sketched a tentative answer, which we will work out in the present section for the 
examples of the scalar and the electromagnetic field. Let us review the basic idea of the procedure 
before we spell out the details: We are going to replace (|5.1| ) by a simpler Hamiltonian which 

• is a good approximation of ( |5.1| ) for slowly varying q and p and 

• is simple enough such that the EOM can be solved exactly. 

The resulting theory will be an approximation for low energies, the detailed information contained 
in the full Hamiltonian (|5.1|) which is only relevant for processes of very high energy gets integrated 



out. This idea underlies also the works [E5| and [26j and, at a rather simple level, is the basis for 



the recovery of continuum elasticity theory from the atomic description in solid state physics (see for 
example p9| ). 



We will now turn to the scalar field Hamiltonian and explain the steps we will take to implement 
the above idea in detail. The Maxwell Hamiltonian will be treated along the same lines in section 



5.2 



5.1 Dispersion Relation for the Scalar Field 

The basic field variables underlying the quantization of the scalar field in 



4> v = <t>{%{ v ))i an d 7T« = / 7T, (5.2) 

J R v 

were represented by the operators — ilnU(v), Y v . R v is the cell containing v in a polyhedral decom- 
position of E dual to 7. According to what we have said in the introduction to this chapter, in the 
considerations to follow we will replace these operators by their classical counterparts ( |5.2j) upon 
assuming the quantum fields to be in a coherent state. 

Using the results of section |4], the Hamiltonian for the scalar field we are considering can be 
written as 

Hkg = 7^— (^kin(vr) + F dcr (0) + K 2 F m ) , (5.3) 
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where 



Pder 



J2 Vdet P 

V 



v/t32 



TrP- 2 (t»] 



!> 2 



+ 



^/det P(u) 



kin 



< Vdet P( 



1 + 



!l73 Tr(p- 2 )P / 2 / ,(^) + ^ 



128 



32 



dt ,(t)vd, 



Col 



e„, T , V«i 



1707 



t 512 



TrP~ 



7T„, 



Now we will express the field quantities 4> v , n v by the basic fields (f>(x), tt(x), using an approximation 
which is good in the case 4>(x),tt(x) vary only very little on the scale e of the graph. The idea is to 
isolate the rough structural properties of ( |5.3|) that lead to corrections as compared to the standard 
dispersion relations and to discard the microscopic details that will only yield higher order corrections 
which are not visible in the long wavelength regime. 

To this end, we Taylor expand the field variables <f> v , tt v around the location x(v) of the vertex v, 
i.e. we make the replacements 

(j) v — ► <f>(x(v)), 

tt v — ► 7r(x(v))Vo\(R v ) + a {a) (v)d a iT(x(v)) + . . . , 

— bfd^X (V)) + bfb^dadaf^v)) + .... 



0. 



<-i 



and truncate the right hand sides at the desired order. Note that in the above formulae we have 
introduced the geometric quantities 



a {a) (v) 



bf{v) = x\f{e I {v)))-x a {v), 
and let us furthermore define 

b { I a) (v) = ^(x a [f(e I (v))]-x a [f(ej(v 

which we will have opportunity to use below. Also, it is perhaps worthwhile to remind the reader at 
this point that all edges are taken to be outgoing from v. 

Then we replace the coefficients of the continuum fields by graph averages and the sums by 
integrals. As argued in |l|, this is a good approximation, as long as <fi and 7r are slowly varying on 
the graph scale e. Let us detail this step for the example of the mass term. We write 



Pm = ^Vol(P„ 



v/detP( 
Vol(P^ 



f 1 

l + -^-TrP- 2 ( 
v^32 V 



5>*Vol(P„ 



^detP(-) f p 1 ,, yggP 
Vol(P.) v^32^ Vol(P. 



TrP- 
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Voi(p.) » + v!32« VoK^o » 

where (( • )) denotes the graph average 

«c(-)» = ^E c w 

for vertex dependent quantities C(v). N denotes the number of vertices of the graph. In case the 
graph has a countably infinite number of vertices, the above definition has to be replaced by an 
appropriate limit of finite sums. 

Analogously, we make the replacements 

-^kin - 



* J (A + Ai) n\x) + A^dM^da^x) + . . . d 3 x, 

F der — > jf + si a)(o,) ) d a <KZ)d*m 



y (Co + C 1 )0 2 (f) + ... cf 3 a: 



where the coefficients in the kinetic term are defined to be 

_ 1707^ V*TrP(v) ^ 
1 512 t KK v /detP(v) /h 

A/det P(-u) 

the ones in the derivative term as 

s w(-') ((Vd^j ty(p- 2 )pa,(.) + &?&?;», 

and finally the coefficients in the mass term by 

C = ((y/detP(v))), 

38 



i el 



Note that we have just written down the leading order terms and the first order corrections, where 
a "first order correction" is either of 

1. A term that is next to leading order in the Taylor expansion and leading order with respect to 
the fluctuation calculation. 

2. A term that is leading order in the Taylor expansion and next to leading order in the fluctuation 
calculation. 

Terms that are leading order in the fluctuation calculation carry a superscript while that are first 
order corrections in the fluctuation calculation are marked by a superscript 1. Finally note that 
we have dropped terms that end up to be a total derivative and therefore do not contribute to the 
Hamiltonian. 

Before we write down the resulting dispersion relation we invoke our restriction to random pro- 
cesses which imply Euclidean invariance on large scales of the resulting random graphs. That is we 



assume 



A 



(a)(a') 



5° 



B 



(a)(a') 



8° 



B (a)(a>) _ 6 a 



For the tensors of fourth rank, the situation is slightly more complicated: 5 ab d~ cd , 5 ac 8 db and S ad S bc 
span the space of rotationally invariant tensors of fourth rank. But contraction of any of them with 
k a k b k c k d is equal to \kf. 

We can now write down the dispersion relation for the Hamiltonian resulting from the above 
replacements 



uo\k) =K 2 [AoCo + Aad + AtCo] 



+ Ifcl 



+i*r 



+ . . . . 



A B + A B^ l) + A X B®V> + K*A^C 



\a b 



(ii)(n) 



(iii)(i) + ^ WjBi 



(i)(i) nW(l) 
o 



(5.4) 



We will discuss the physical content of (5.4) in section |5.3| . Before that, we give a similar calculation 
for the electromagnetic field. 



5.2 The Electromagnetic Field 

This section is devoted to the calculation of an (approximate) dispersion relation for the electro- 
magnetic field. The treatment is completely analogous to the one given for the scalar field in the 
last section, so we can be rather brief, here. Again, we introduce the continuum fields A a (x), E a (x) 
underlying the regularization and quantization of the Hamiltonian performed in || and also the 
classical quantities 



E e = E which was represented by Y e , 

Jse (5.5) 
A e = I A which was represented by — i In H e 
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Figure 1: The loop ai(v) 



(subject to the subtleties associated with the logarithm spelled out in detail in in the quantum 
Hamiltonian (|4.50 ). When we replace the gravitational operators in the Hamiltonian by their expec- 
tation values obtained in the last section and the operators for the matter fields by their classical 
counterparts (|5.5f) we get 



{Hm,~{) 



2Q 



(P el (P) + P mag (P)) 



where 



M 



v Ial'a' 



v IV 



t V512 

4EE [ V^PRP- + f (jg^ 2 Tr P- - Hp- 



16 



^E e ^(v)E e ^(v), 



A~ A~ 

■"-Oil f^-CLjl 1 



where a(v) is the loop around the vertex v "in the /-plane" as depicted in figure [I]. Now we Taylor- 
expand the A a , E e . To this end we introduce some geometric quantities: 



s e Jv) 



n a {y) dy, 



n a {y){y- x(v)) b dy, 



n a {y){y-x{v)) b (y- x(v)) c dy, 



where n denotes the normal to the surface of integration. Moreover 



sl(v) 



sVM - Sc 1 {V, 



sib(v) = 2 ( S ab( W ) - S lb ( V ) 



Now we can make the replacement for the electric field: 

E e > s e a (v)E a (x(v)) + s e ab (v)d b E a (x(v)) + .... 

We proceed in a similar fashion for the connection: 



a a (s)(a(s) - x{v))\ b a a bc (v) = / a a (s)(a(s) - x(v))\a(s) - x(v)) 



(5.6) 
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whence we replace 



A a ► b ab (v)d b A a {x(v)) + ^bf c (v)d b d c A a (x(v)) + .... 



Inserting this into the expressions for F el and F mag and subsequently replacing the resulting coeffi- 
cients by graph averages results in the total replacement 

F* > I {S$w + 5 S(a ) E a (x)E a \x) + 2S^ )(albl) E\x)d b ' E«\x) + ..., 

F mag > (B^ a ' hl) + B ^ a ' b 'A d b A a {x)d v A al {x) + 5S a6 ^ a ' 6 ' c '^ 6 A a (f)a 6 ^ c ,A ,(^) + • • • , 



where 



77' 
?4 



77' 

^(a)(a'6') = E ((>^^( • " " )», 

77' 

and analogously 

5r )(aV) =E<<^r,(-)^(-)C(-)>>, 

77' 11 

! 4 PV ^,/ 1 /763 „ _„ N 13 



B WM = ^((y {^7? Tr P-\ ■ ) - ^PfA ■ • • ))», 

77' u 



^ W = E «^^( • ■ K d { • )»• 



77' 



Now we can make the replacements ( |5.7| ) and obtain a Hamiltonian for the continuum fields A a (x) 
E a (x). A straightforward calculation yields the resulting equations of motion: 



A d =d h d h ,A a 'S [d){a) B^t^) - B^ db '^ 

+ <9Ad c >' [flWC'O - W) + (5 (d)(ac0 - S {a)(dd) ) [B^ a ' b ^ - 
+ ... 

M , oC 1 ) 



(5i 



where we have used shorthands 5 , ( a )(a') = ^( a j( a ')+^(a)(a') an< ^ srmnar ly f° r the other tensors Sqq, 

Before we spell out the resulting dispersion relation, we use the rotation invariance of the graph 
on large scales: It is clear that 

S(a)(a') ~ Sao,', S(a)(a'b') ~ ^aa'b'- 

For the tensors of higher rank, the situation is slightly more complicated: For rank four, the space 
of invariant tensors is three dimensional, the space of rank five tensors is ten dimensional. But if we 
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take into consideration the symmetries of the terms, these tensors get contracted with, there is only 
one invariant tensor left in each case. We define 



(0/1) j_ 1 qCV 1 ) „(0/l) j_ 1 ( p>(ab)(ab) _ ST^ R (aa)(aa) \ 

C l - 3 Z^ D (a)(a)' C 3 "g Z^V 2^^0/1 I' 

i \ ab a / 

„ j_ 1 abc q(P) „ ^ 1 \^ f: Jj{bc){acc) 

C 2 - g 2^ 6 (a) (be)' C 5 - g efeaCjW 

abc abc 

A straightforward calculation shows that the equations of motion (|5.8|) simplify to 

2(t, s) = (ci 0) cf + cf } 4 1} + c[ 1} 4 0) ) AA(t, x) + (c 2 cf - c[ 0) c 5 ) A rot A(t, x). (5.9) 

Note that in the last equation we have just kept terms of leading order and first order corrections, in 
the sense that we have explained in the previous section. Also we have eliminated a term containing 
div A by choosing the appropriate gauge. 

Equation ( |5.9| ) leads to a chiral modification of the dispersion relation for electromagnetic waves: 
Let a unit vector e 3 be given and choose ei,e*2 such that the form a righthanded orthonormal 
triple. Then a circularly polarized wave of helicity ±, propagating in the direction given by e% can 
be written as 

A k (t, x) = A [e*i cos {u)±{k)t — k e 3 ■ x) ± e 2 sin {u)±{k)t — k e 3 • x)] . 
This is a solution to the wave equation ( p.9| ) provided that 

co ± (k) = \k\ ^(cfcf + ci 0) 4 1} + c^cf) ± (c 2 4 0) - cS 0) c 5 ) k. (5.10) 

Thus we have found a chiral modification to the dispersion relation. Note that this chiral modification 
is similar but not completely analogous to the birefringence occurring for light propagation in some 
crystals. The latter effect is not isotropic, it also depends on the direction of propagation relative to 
the symmetry axes of the crystal, whereas the chiral effect found here is isotropic. This can be seen 
from the fact that nothing in the above formulae depends on the direction of the vector e*3. We can 
now proceed to a discussion of results. 

5.3 Discussion 

Let us start the discussion of the results of the last section by considering the physical units and 
orders of magnitude of the various terms appearing. We will use Fd er , the derivative term in the 
scalar field Hamiltonian, as an example - similar considerations apply to the other terms. 

The classical term corresponding to B^^ ^ + B[ a ^ a ^ is a/ det qq aa ' . The latter is dimensionless, 
since q is. B^^ a ■* has the structure 

~ -L-fL=bb (5.11) 
Volv/detP 
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where Vol is a volume. Since [P] = meter 2 , [P 2 /-y/det P] = meter, b is also a length, unit-wise, so 
£W( a ) j s indeed dimensionless. B[ a ^ a ^ has the structure 

*<■><■> -A^Zl-y (5.i2) 

1 tVol VdetP 

so it is again dimensionless as it should be. The structure of B^ h ^ a h ^ is 

~ ^^L=bbbb, (5.13) 
Volv/detP 

so its unit is meter 2 which is the correct one for a term proportional to \k\ 4 in the dispersion relation. 
As for orders of magnitude, we remark the following. Assume q a b = 0(1) in the chosen coordinate 
system. Then 

p = 0(e 2 ), Vol = 0(e 3 ) and b = 0(e). (5.14) 
Using ( |5.11| ) it follows that B^ a ^ a ^ = 0(1), so the leading order term has the right order of magnitude. 



As for the order of magnitude of B[ a ^ a \ we use (|5.12|) and ( |5.14| ) to conclude that 



t e 




which is very small since a < 1/2. 

Consider finally B^ b)(a ' bl) : From (|[TJ) and (^141) we see that B^ b){a ' b ' ] = 0(e 2 ) 



As for the other terms in the dispersion relation, similar results can be seen to hold: The leading 
order term has same unit and order of magnitude as the corresponding classical term and the ratio 
of leading order to first order correction is of order t 1 ~ 2a '. 

We will now discuss the structure of the dispersion relations ( |5.4p , ( [5.101 ). The coefficients ap- 



pearing are given as graph averages of certain local geometric quantities of the random graph. Let us 
call these graph averages moments of the random graph prescription (RGP for short). So, in order 
to get numerical statements from the results of the last section, one has to fix the scale L, an RGP 
and compute the relevant moments. Such a computation might be hard to perform analytically, 
but a computer can easily determine the moments occurring in ( |5.4| ),( [5"Tl~0| ) for a given RGP, so this 
calculation does not present a principal difficulty. 

The more serious issue here is that there are certainly many RGPs, all leading to different graph 
averages and hence different predictions, and it is a priori not clear how one can single out the 
"right" one. We note however that for not too pathological RGPs, the graph averages should be 
approximately equal so that at least the size of the different terms in the dispersion relations is not 
too sensitive on the choice of the RGP. Moreover, again for a not too pathological RGP, the moments 
showing up in the dispersion relations should be related. To give an example, a plausible assumption 
is that 



«^K«vsb») 
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and that their difference would not depend very strongly on the chosen prescription. Thus there will 
be approximate relations between the different coefficients in the dispersion relations which are not 
affected by the choice of a specific RGP. 

Moreover we note that the leading order terms in the coefficients depend on the RGP. This might 
at first seem to be a problem as well, since it means that we will have to tune the RGP in such a 
way that the leading order terms assume their classical values. On the other hand, this can be seen 
as a blessing: Fixing the leading order term means to fix one moment of the RGP. Via the relations 
conjectured above, this will also approximately fix other moments, independently of the specific RGP 
assumed, and thereby to a certain extent the higher order corrections. 

Investigations in this direction are worthwhile but beyond the scope of the present work. Let us 
for the rest of this section assume that a prescription is fixed and the relevant graph averages have 
been computed. 

Next we observe that two different sorts of corrections appear in the dispersion relations: The 
first sort of correction is simply a correction to the leading order term coming out of the fluctuation 
calculation of section Its relative magnitude was found to be t x ~ 2a . We will call this sort of 
correction a fluctuation correction. 

The other sort of correction is a term containing a higher power of \k\ as compared to the standard 
dispersion relation. We will call this kind of correction a lattice correction. We have demonstrated 
for the example of BQ^ ab ' that the terms proportional to |/c| 4 are of the order e 2 , therefore the 
relative magnitude of the lattice corrections is of the order 

°(£) = ^ °(n<o(n 

3 

Similarly the terms proportional to \k\ in the dispersion relation for the electromagnetic field are of 
the order t a L/X. 



When comparing our results for the electromagnetic field with the ones of [25. 26 1 we find the 



following: The result of Pullin and Gambini |25j does not contain any fluctuation corrections. This 
is however not result of the calculation but rather assumed from the beginning. As for the lattice 
corrections, they find a chiral modification to the dispersion relation as we do here. The relative 
magnitude of the correction is however £p/X. 



Alfaro et al. [26] also do not have fluctuation corrections by assumption. They find the helicity 



dependent correction of [25| and the present work, again of the order £p/X. They also get higher 
order corrections the precise structure of which depends on a parameter which is not fixed. 



Thus our results agree with that of |p5|, |26| as far as the structure of the dispersion relation 
is concerned. We additionally have fluctuation corrections and, most importantly, the corrections 
found do not scale with an integer power of £p, contrary to their finding. This signals a warning to 
assumptions made in |^0[ to take into account only corrections which are of the order (£p/L) n where 
n is an integer. Notice that the fluctuation correction and the lattice correction are equal at a = ~. 
Thus the leading correction is always of the order of at least t 1 / 3 . 

Finally we should make a few remarks concerning a possible detection of the corrections in 
experiments. The fluctuation corrections will not show up in an experiment testing for a frequency 
dependence of the velocity c of light, since they merely correspond to a frequency independent shift 
of c. Also, these corrections are certainly not measurable by measuring the flight-time of photons 
since their velocity would already be the "bare" leading order term plus the fluctuation correction. 
Fluctuation corrections may however be measurable by comparing flight-times of photons in different 
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geometries, since the corrections will change when the calculations presented in this chapter are 
repeated with LQC approximating a non-flat spacetime. To discuss how this could be done in 
practice is however beyond the scope of the present work. 

Whether the lattice corrections are big enough to be detectable in the data from current or planned 
7-ray burst observations crucially depends on the values of a and L. For the value a = 1/3 which 
renders fluctuation and lattice corrections equal in magnitude (and which is close to the lower bound 
value 2/5 derived in 0), and L of the order of a 7-ray wavelength, a rough estimate shows that the 
lattice corrections would indeed be detectable in the foreseeable future. 

So, to conclude this chapter, we should repeat that not too deep a significance should be attached 
to the precise values of the coefficients in the dispersion relations obtained: There are still some 
ambiguities present in the GCS which we will discuss below: The quantization of the Hamiltonians, 
in the procedure to obtain the dispersion relations from the expectation values and, as a consequence, 
in the coefficients themselves. Also the replacement SU(2) — > U(l) will certainly affect the precise 
numerical outcome. Most significantly, so far we have little control on what will happen to the size 
of quantum corrections when our kinematical coherent states are replaced by physical ones. Within 
our kinematical scheme the structure of the corrections, as well as the orders of magnitude t a , t l ~ 2a 
of the two sorts of corrections are robust, however. Thus we are possibly in trouble because such 
corrections seem to lie in the detectable regime. If such corrections are not found, then presumably 
it is not justified to use kinematical coherent states. 

Finally, the approximate relations between the different graph averages will make the predictions 
of a more complete calculation much less dependent on the random graph prescription chosen, then 
one might at first fear. Similar remarks apply if, as advocated for example in ||, instead of working 
with a fixed random graph, one averages over many of them. (Notice that also in that case, averaging 
procedures are not unique). In order to remove those ambiguities one should probably set up a 
variational principle in order to optimize a family of semiclassical states according to a given set of 
observables. 



6 Summary and Outlook 

In the this work we have presented a calculation of dispersion relations for the scalar and the elec- 
tromagnetic field coupled to quantum general relativity. These dispersion relations bear corrections 
to the standard ones, due to the discreteness of the states of the geometry and to the bound on the 
uncertainty product of configuration and momentum variables in QGR. The calculations rest on the 
quantization of the matter parts of the Hamilton constraint given in |l| and the coherent states for 
QGR constructed and analyzed in || |H], [□], [19j 



Corrections to dispersion relations due to QGR were also computed in [^, g6j and the present 
work partly rests on the ideas implicit and explicit in these pioneering works. The form of the 
correction term in the dispersion relation for the electromagnetic field found in the present work 
agrees with that of PH |2l)fl . This is not too big a surprise since there is no other rotation invariant 
term in k of the same order. However, we find important differences in the order of magnitude of the 
effects, as compared to ]23], |2l)fl . Moreover, the results of the present work are more specific, since a 
definite class of semiclassical states, the coherent states for QGR are employed in the calculation. 

Rather than making precise numerical predictions, the aim of the present work is to demonstrate 
the steps necessary in such a calculation, to highlight the issues that remain to be clarified and to 
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give a robust estimate of the size of the effects. 

In this spirit, we have simplified the calculation of the expectation values in ^ by replacing the 
full gauge group by its Ionii-Wigner limit U(l) 3 . This replacement will certainly affect the precise 
numerical outcome but not the order of magnitude of the correction. Also we have not specified a 
prescription for obtaining random graphs, but only assumed general properties that such a procedure 
will have. Most importantly, the effect of using kinematical rather than physical coherent states is 
presently not well understood. 

The main achievements of the present work can be summarized as follows: 

The calculation given in section [| shows how expectation values of complicated operators in 
coherent states for quantum general relativity can be computed and there is no principal difficulty 
in repeating such a calculation for the full gauge group SU(2). 

Perhaps even more important are the order of magnitude estimates of the resulting effects obtained 
in this work: They depend on very few parameters and will continue to hold true when more general 
complexifier coherent states are used. The main choices that enter are: 

• A complexifier C has to be chosen for the construction of the coherent states. (Of course, there 
are more general semiclassical states than coherent ones). 

• A class of observables has to be chosen that should be approximated well by the coherent states. 

• A (random, averaged) graph has to be chosen. 

The other parameters are fixed by the above choices: The requirement that C/his dimensionless 
forces the parameter t in the definition of the resulting coherent states to be (£p/a) n where n is some 
positive number and a a length scale which is not yet fixed. 

The nature of these observables (do they involve one-, two- or three-dimensional integrations? 
etc.) determines a) a length scale L and b) the exponent j3 in the expression for the classical error 
(e/Lr. 

The length scale a gets fixed to be L by requiring fluctuations of configuration and momentum 
degrees of freedom to be equal. Finally, the typical edge-length e of the random graph is found to 
be a weighted geometric mean by requiring the fluctuations to be minimal. Thus, at least within the 
vast class of complexifier coherent states, the structure of the ambiguities and their principal effects 
on the orders of the magnitudes of the quantum corrections can be neatly classified] 

Many things remain to be done before one can really obtain reliable predictions of observable 
effects from quantum general relativity: 

The procedure used to obtain dispersion relations from the discrete classical Hamiltonians has to 
be further analyzed, and rigorously justified at least in models which can be solved analytically. The 
influence of the choice of a random graph should be investigated, and concrete procedures have to 
be implemented. A more distant goal is to also analyze possible back reaction effects of the matter 
on the gravitational field. These were neglected in M and in this work since it would require to solve 
the combined matter - geometry Hamiltonian constraint and force us to work with physical coherent 
states. 

Thus, although we certainly did not carry out a first principle calculation, we hope to have made a 
modest contribution to an understanding what the principal problems are and how such a computa- 
tion could possibly be carried out in principle. Also, we hope to have demonstrated that QGR is still 
far from making reliable semiclassical predictions until one is convinced of the physical relevance of a 
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definite scheme. However, it should have become clear that once such a scheme has been identified, 
QGR is able to provide precise numerical predictions. In any case, at least for the limited purpose of 
showing that some version of the quantum Hamiltonian constraint is correct (for which kinematical 
coherent states are unavoidable), the results of the present two papers should be relevant. 
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